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PREFACE 

The  Rand/ARPA  Climate  Dynamics  Project,  sponsored  by  the  Defense 
Advanced  Research  Projects  Agency,  uses  computer  simulations  of  the 
earth's  climate  as  a  means  of  understanding  and  predicting  climate  and 
its  sensitivity  to  modification  by  human  action,  whether  deliberate  or 
inadvertent.  The  basic  research  tool  is  a  Rand-adapted  version  of  the 
Mintz-\rakawa  (M/A)  two-level  atmospheric  general  circulation  model  de¬ 
veloped  at  the  University  of  California,  Los  Angeles,  under  National 
Science  Foundation  sponsorship  (see  A  Documentation  of  the  Mintz- 
Arakaua  Two-Level  Atmospheric  General  Circulation  Model ,  Rand  Report 
r_877-ARTA,  by  W.  L.  Gates,  E.  S.  Batten,  A.  B.  Kahle,  and  A.  B.  Nelson). 

A  characteristic  of  the  M/A  simulations  is  the  production  of  large 
volumes  of  data.  One  60-day  simulation  now  produces  roughly  fifteen 
million  numbers.  This  report  documents  several  computer  programs  that 
handle  this  massive  output,  reducing  it  to  more  manageable  and  compre¬ 
hensible  form:  three  that  generate  maps,  profiles,  and  tables;  one 
that  calculates  and  plots  certain  differences;  and  two  auxiliarv  pro¬ 
grams  designed  for  more  convenient  and  efficient  use  of  the  main  pro 
gram. 

The  author  would  like  to  thank  Robert  Mobley  and  E.  S.  Batten  for 
reviewing  an  earlier  draft  of  this  report,  and  Marian  Winston  for  her 
editorial  help. 
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SUMMARY 


This  report  is  a  user's  guide  to  six  computer  programs  employed  by 
the  Rand/ARPA  Climate  Dynamics  Project  in  connection  with  simulations 
of  the  global  atmosphere  to  study  the  highly  complex  interactions  that 
produce  world  climate.  The  primary  research  tool  of  the  project,  the 
Mintz-Arakawa  two-level  general  circulation  model,  is  a  massive  program 
requiring  25  hours  actual  CPU  computation  time  for  a  60-day  simulation 
and  generating  two  billion  numbers,  of  which  15  million  are  stored  in 
18  large  arrays. 

The  MAPGEN,  Map  Display,  and  Zonal  Profile  programs,  documented 
here,  do  postprocessing  of  this  output,  manipulating  it  into  more  man¬ 
ageable  and  comprehensible  form — maps,  profiles,  and  tables.  The  Dif¬ 
ference  program  calculates  the  differences  between  experiments  in  a 
form  that  car.  be  plotted  by  the  Map  Display  program.  Two  auxiliary 
programs  are  documented  here  for  completeness:  Quick  Look  and  Disk 
Initializer,  which  greatly  facilitate  use  of  the  M/A  program. 

The  MAPGEN  and  Map  Display  programs  are  run  one  after  another  to 
produce  contour  maps  of  the  atmospheric  variables.  Together  with  the 
Zonal  Profile  program,  they  create  a  standard  set  of  displays  of  the 
atmosphere  over  the  earth's  land  masses  and  oceans.  MAPGEN,  a  very 
large  program,  is  used  to  obtain  values  and  averages  of  a  set  of  mete¬ 
orological  variables  that  are  functions  of  the  primary  dependent  vari¬ 
ables;  up  to  28  variables  can  be  obtained,  some  of  them  at  both  of  two 
atmospheric  levels,  for  each  of  the  3,312  nodes  on  the  46  *  72  latitude/ 
longitude  grid.  These  values  are  put  on  tape  for  input  to  Map  Display. 
MAPGEN  consists  of  about  4000  FORTRAN  IV  source  cards,  uses  1,250,000 
bytes  of  core,  and  requires  50  minutes  to  process  the  output  from  a 
60-day  simulation. 

Map  Display  uses  about  150,000  bytes  of  core  and  about  20  seconds 
of  computing  time  per  contour  plot.  A  contour  map  and  a  profile  map 
of  the  zonal  means  are  output  for  each  set  of  values  generated  in  MAP¬ 
GEN.  The  standard  output  package  consists  of  50  different  contour  maps 
and  the  zonal  mean  plots  for  each  variable,  accompanied  by  tables  of 
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the  values ,  the  zonal  means,  the  meridional  means,  the  zonal  means  over 
oceans,  and  the  global  mean.  Very  little  computing  time  is  required 
for  the  profile  plotting  routines. 

Since  the  standard  M/A  output  uses  six-hour  accumulations,  and  cer¬ 
tain  variables  change  too  rapidly  to  be  satisfactorily  represented  by 
six-hour  samplings,  the  Zonal  Profile  program  works  from  the  accumula¬ 
tion  of  half-hour  values  for  the  radiation,  evaporation,  and  heating 
components.  This  program's  output  package  consists  of  up  to  32  profile 
graphs.  Zonal  Profile  consists  of  about  1500  source  cards  and  uses 
about  150,000  bytes  of  core  and  18  seconds  computing  time  for  a  60-day 
simulation. 

The  Disk  Initializer,  Difference,  and  Quick  Look  programs  are  not 
formal  output  programs,  but  were  written  to  facilitate  the  running  of 
'.he  M/A  model.  These  are  much  smaller  programs.  Disk  Initializer  sets 
up  a  new  disk  for  either  a  new  or  continuation  experiment.  Quick  Look 
can  be  used  to  stop  the  M/A  simulation,  read  the  disk,  and  print  out 
the  values,  to  enable  the  researcher  to  determine  if  the  M/A  program 
is  operating  as  expected— a  particularly  useful  function  with  a  program 
that  runs  so  many  hours.  The  Difference  program,  as  mentioned  above, 
expresses  the  difference  between  two  simulation  experiments  in  a  form 
that  can  be  plotted  by  Map  Display. 
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I .  INTRODUCTION 


This  report  is  a  user's  guide  to  six  computer  program.:  used  in  the 
Rand/ARPA  Climate  Dynamics  Project.  Under  the  sponsorship  of  the  Ad¬ 
vanced  Research  Projects  Agency,  the  project  performs  computer  simula¬ 
tions  of  the  global  atmosphere  to  study  the  highly  complex  interactions 
that  produce  world  climate.  Human  activities  may  have  already  produced 
significant  changes  in  the  atmosphere,  yet  our  knowledge  of  the  dynam¬ 
ics  of  climate  is  not  advanced  enough  to  permit  us  to  predict  with  con¬ 
fidence  the  results  of  such  changes. 

The  primary  research  tool  of  the  Climate  Dynamics  Project  is  the 
computer  simulation  of  world-wide  circulation  patterns,  using  a  version 
of  the  Mintz-Arakawa  (M/A)  two-level  general  circulation  model  origin¬ 
ally  developed  at  the  University  of  California.  The  basic  Rand  adapta¬ 
tion  of  the  M/A  model  is  documented  in  Ref.  1,  which  should  be  consulted 
f jr  more  complete  descriptions  of  the  physical  variables  and  the  equa¬ 
tions  used  in  the  simulation.  There  have  been  changes  made  in  the 
model  output  since  that  report  was  published,  and  some  of  these  are 
documented  here. 

The  output  now  comprises  18  arrays  storing  a  total  of  15  million 
numbers  per  60-day  simulation,  rather  than  the  10  arravs  stoied  pre 
viously.  Three  new  programs — MAPGEN,  Map  Display,  and  Zonal  l-'ofile 
do  postprocessing  of  the  M/A  output.  Three  auxiliary  programs— Quick 
Look,  Difference,  and  Disk  Initializer — facilitate  the  use  of  the  M/A 
program.  An  older  version  of  MAPGF.N,  descril  .1  in  Ref.  1,  is  super¬ 
seded  by  the  version  described  below. 

THE  MJNTZ-ARA^V^  PROGRAM 

The  basic  numerical  formulation  concerns  the  division  of  the 
earth's  atmosphere  into  a  grid  defined  by  \  degrees  of  latitude  and  ) 
degrees  of  longitude.  Thus,  there  is  a  W  '  72  array  with  1,112  nodal 
points.  For  each  node,  the  state  variables  of  the  atmosphere  are  de- 
rined  at  each  of  two  levels  defined  by  a  sigma  surface  which  is  propor- 

The  values  saved  are  calculated  at  the  sigma 


tional  to  pressure. 


surfaces  of  one-quarter  and  three-quarters,  where  sigma  equals  zero  at 
the  top  of  the  troposphere  (200  mb)  and  1  on  the  earth's  surface.  Se¬ 
lected  meteorological  variables  (discussed  more  fully  below)  are  com¬ 
puted  for  each  node  point  every  simulated  half-hour.  At  six-hour  in¬ 
tervals,  the  state  variables  at  that  time  point  are  computed  and  stored 
in  the  arrays  mentioned  above. 

The  simulation  of  a  60-day  period  generates  approximately  2  billion 

numbers,  of  which  15  million  are  stored  on  disk,  and  requires  25  hours 

CPU  computing  time  on  the  IBM  360  91,  one  of  the  fastest  computers  in 
* 

existence.  The  MAPGEN,  Map  Display,  and  Zonal  Profile  programs  manip¬ 
ulate  this  massive  output,  reducing  it  to  manageable  and  comprehensible 
form,  and  the  Difference  program  can  be  used  to  plot  the  comparison  be¬ 
tween  two  simulation  experiments. 

THE  M/A  OUTPUT  PROGRAMS 

Certain  output  routines  were  included  in  the  original  Mintz- 
Arakawa  two-level  general  circulation  model  as  obtained  from  UCLA.  Us¬ 
ing  an  overlay  technique,  the  program  accumulated  the  output  arrays  as 
it  ran,  and  after  30  simulated  days,  the  output  arrays  were  written  and 
a  contour  map  was  printed.  At  Rand,  A.  B.  Nelson  deleted  the  original 
output  option  and  wrote  the  first  versions  of  MAPGEN  and  Map  Display. 
The  present  versions  and  the  Zonal  Profile  program  were  written  by  the 
author,  under  the  direction  of  E,  S.  Batten  and  A.  B.  Kahle.  The  ver¬ 
sions  described  here  have  been  in  use  since  February  1972. 

MAPCEN  is  used  to  obtain  values  and  averages  of  a  set  of  meteoro¬ 
logical  variables,  such  as  radiation,  cloudiness,  rain,  evaporation, 
albedo,  etc.,  which  are  functions  of  the  primary  dependent  variables. 
These  values  are  then  put  onto  tape  for  the  Map  Display  program.  In 
the  present  version,  up  to  28  different  variables  can  be  obtained,  with 
some  of  these  available  at  more  than  one  level.  The  standard  output 
package  consists  of  50  different  maps.  Map  Display  is  used  to  generate 
contour  maps  from  the  MAPGEN  output.  Using  IGS  S-C  4060  graphics 

it 

A  modified  version  of  the  program  now  requires  15  hours  CPU  time. 
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routines  (Ref.  2),  contours  of  the  data  are  drawn  over  a  background  map 
o f  the  world,  and  the  plots  of  zonal  means  are  also  drawn  for  each  var¬ 
iable.  The  values  themselves  are  also  printed  in  tabular  form,  together 
with  the  zonal  means,  the  meridional  means,  the  zonal  means  over  oceans, 
and  the  global  mean. 

Since  MAI  .N  samples  computed  variables  once  every  six  simulated 
hours,  the  normal  daily  variation  of  some  of  these  variables,  such  as 
radiation,  evaporation,  and  heating  components,  can  cause  a  significant 
sampling  error  in  the  results.  For  studies  of  heat,  radiation,  and 
moisture  balance,  it  is  desirable  to  eliminate  this  error.  Therefore, 
to  obtain  more  accurate  zonal  means,  the  present  version  of  the  M/A 
program  accumulates  zonal  arrays  each  half-hour  for  these  variables  and 
stores  the  values  on  disk.  The  Zonal  Profile  program  was  written  to 
collect  and  plot  zonal  profiles  of  these  variables. 

The  MAPGEN  Program 

MAP GEN  is  a  very  large  output  program.  The  program  consists  of 
about  4000  FORTRAN  IV  source  cards,  and  uses  1,250,000  bytes  of  core. 

To  generate  the  values  for  the  50  standard  maps  for  a  30-day  accumula¬ 
tion,  MAPGEN  uses  up  to  25  minutes  of  CPU  time,  costing  about  $250  on 
the  IBM  360/91  at  UCLA. 

The  principal  output  from  MAPGEN  is  a  tape  of  the  time-averaged 
values  and/or  snapshot  values  (arrays  at  one  time  point) .  Twenty-eight 
different  variables  can  be  generated. 

The  Map  Display  Program 

The  Map  Display  program  uses  approximately  150,000  bytes  of  core 
and  about  20  seconds  of  CPU  time  per  contour  plot  on  the  360/91.  The 
printing  and  profile  plots  use  very  little  CPU  time. 

For  each  set  of  values  generated  in  MAPGEN,  a  contour  map  and  a 
profile  plot  of  the  zonal  means  are  generated  in  Map  Display.  The  ar¬ 
ray  values  are  also  printed,  together  with  the  zonal  means,  the  zonal 
means  over  oceans,  the  global  means,  the  meridional  means,  and  the  map 
title.  The  contour  map,  identifying  the  local  high  and  low  values  of 
the  array,  is  drawn  on  a  base  map  outlining  the  continents. 


The  values  of  the  contours  drawr  are  determined  by  specifying  the 
contour  Interval  and  the  value  of  one  contour  to  be  drawn  as  a  dashed 
line.  The  otiier  isoline  values  are  positive  and  negative  multiples  of 
the  isoline  Interval  from  the  dashed-line  value.  Alternatively,  the 
user  may  input  up  to  15  arbitrary  isoline  values,  including  his  choice 
of  the  dashed-line  value.  See  Fig.  1,  p.  29,  for  a  sample  contour  map. 

The  contour  map  is  usually  a  rectangular  projection,  but  polar 
projections  can  also  be  plotted  if  desired.  The  polar  projections  can 
be  either  stereographic  or  equidistant  projections,  and  of  course  can 
be  centered  on  either  the  North  or  the  South  Pole. 

For  the  profile  plot  of  the  zonal  means,  the  program  finds  the 
largest  value  of  the  zonal  means  and  uses  that  as  the  scale  factor  for 
the  graph.  While  this  method  of  scaling  has  advantages,  it  makes  com¬ 
parison  of  these  graphs  difficult.  (See  Fig.  2  on  p.  30  for  a  sample 
zonal  means  graph.) 

Map  Display  output  options  enable  the  user  to  select  the  printout, 
contour  map,  and/or  zonal  means  graph. 

The  Zonal  Profile  Program 

The  Zonal  Profile  Program  uses  approximately  150,000  bytes  of  core 
and  requires  nine  seconds  of  CPU  time  to  generate  32  profile  graphs 
from  a  simulated  30-day  period.  The  program  consists  of  about  1500 
FORTRAN  IV  source  cards.  A  sample  profile  graph  is  given  in  Fig.  3  on 
p.  34 . 

From  th«  zonal  mean*-  saved  during  the  integration  of  the  M/A  pro¬ 
gram,  the  Zonal  Profile  Program  generates  up  to  32  different  profiles, 
including  profiles  of  radiation  components,  evaporation,  cloudiness, 
and  energy  balance.  The  zonal  variables  are  averaged  over  a  given  f*me 
period  from  values  calculated  each  half-hour,  and  then  are  output  in 
both  printed  and  graphic  forms. 

The  Smaller  Programs 

In  addition  to  the  three  formal  output  programs,  three  smaller 
programs  are  necessary  to  use  the  M/A  program  effectively.  These  are 
called  the  Quick  Look  Program,  the  Difference  Program,  and  the  Disk 
Initializer  Program. 


The  Quick  Look  Program  reads  the  M/A  disk  output  and  prints  out  the 
data.  At  checkpoints  throughout  the  M/A  simulation,  the  M/A  program 
may  be  stopped  and  the  data  on  disk  inspected,  via  Quick  Look,  to  de¬ 
termine  if  the  program  is  behaving  as  expected.  Both  the  state  vari¬ 
ables  and  the  accumulated  zonal  means  are  printed  by  the  program. 

The  Difference  program  is  used  to  calculate  the  differences  between 
two  experiments.  Two  data  tapes  generated  by  MAPGEN  are  used  to  form 
a  new  tape  to  be  used  by  the  Map  Display  program  to  tou..  m«p«  of  the 
differences. 

The  Disk  Initializer  program  is  used  exclusively  by  the  programmer 
in  conjunction  with  the  running  of  the  M/A  program.  This  program  will 
set  up  a  new  disk,  either  for  the  continuation  of  an  experiment  or  for 
the  beginning  of  a  new  experiment. 

ORGANIZATION  OF  THIS  REPORT 

The  next  section  describes  the  new  format  of  the  M/A  disk  output. 
This  is  followed  by  user's  guides  for  the  six  programs. 

Four  appendices  are  given  at  the  end  of  the  program.  Appendix  A 
is  a  lifting  of  the  standard  input  package  for  MAPGEN  and  Map  Display. 
Appendic  B  is  a  printout  of  the  data  on  a  simple  variabl?  (degrees 
Centigrade)  from  the  Map  Display  program.  Appendix  C  gives  the  cor¬ 
respondence  between  the  new  and  old  numbers  for  the  MAPGEN  map-types, 
as  documented  here  and  as  documented  previously  in  Ref.  1.  And  Appen¬ 
dix  D  presents  tabular  detail  concerning  program  variables,  symbols, 
and  arrays. 


II.  M/A  DISK  OUTPUT 


To  understand  any  of  the  M/A  output  programs,  the  user  must  first  under¬ 
stand  what  the  M/A  program  has  written  onto  disk  during  its  simulation. 

The  M/A  simulation  program  forms  two  data  sets  for  a  time  history 
of  the  run.  The  larger  of  the  data  sets  consists  of  large  46  *  72  grid- 
point  arrays,  with  the  smaller  consisting  of  a  number  of  46  *  1  zonal 
arrays . 


THE  LARGE  ARRAY  DATA  SET 

At  the  beginning  of  the  large  array  data  set,  two  arrays  that  never 
change  throughout  the  simulation  are  given.  These  are  TOPCG,  the  sur¬ 
face  topography  indicator,  and  ALBEDO,  the  surface  albedo,  which  appears 
as  ALS  in  the  COMP 3  subroutine. 


Table  1 

SYMBOLS  USED  IN  THE  TOPOGRAPHY  ARRAY 


1.  If  TOPOG  >  1.0,  then  the  grid  point 
is  ocean  and  the  value  is  the  ocean 
temperature  in  degrees  Kelvin. 

2.  If  TOPOG  <  1.0  and  >  -106,  then  the 
grid  point  is  land  and  the  value  is 
the  negative  of  the  altitude  in  meters 
multiplied  by  the  gravitational  con¬ 
stant,  GRAV  (9.81). 

3.  If  TOPOG  <  -1.0  x  106,  then  the  grid 
point  is  land  ice  and  the  value  is 
-1.0  x  106  minus  the  altitude  in 
meters  multiplied  by  the  gravitational 
constant. 

4.  If  TOPOG  -  -1.0  x  106,  then  the  grid 
point  is  sea  ice. 


Once  every  six  hours  of  simulated  time,  the  program  writes  data 
onto  the  disk  in  the  two  data  sets.  The  large  array  data  set  can  be 
divided  into  four  sections:  (1)  time  record,  (2)  the  pre-COMP3  values 
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(explained  below),  (3)  the  state  variabler ,  and  (A)  the  precipitation 
arrays.  Table  2  gives  the  arrangement  on  disk  for  each  time  point  for 
this  data  set.  The  other  data  sets  consist  of  zonal  arrays. 


Table  2 

FORMAT  FOR  EACH  COMPLETE  SET  OF  TIME-DEPENDENT  DATA 


TAU, C (800) 


Pre-C0MP3  Arrays 


T (46, 72, 2) 
Q3(46,72) 
GW(46, 72) 
GT(46, 72) 


State  Variable  Arrays 


P(46,72) 

U(A6,72,2) 

V(A6, 72, 2) 

T  (46,72 , 2) 

Q’>(46,72) 

GW(46,72) 

GT(46,72) 


Precipitation  Arrays 


ACPPLS(46,72) 

ACPPC1(46,72) 

ACPPC3(46,72) 


The  Time  Record  (TAU)  and  Constant  Array 

The  first  record  at  each  time  point  contains  the  reference  time, 
TAU,  in  hours,  and  an  array,  C,  of  800  constants  used  in  the  M/A  simu¬ 
lation  and  the  output  programs. 

For  historical  reasons,  the  starting  point  of  the  Rand-M/A  simula¬ 
tions  is  not  TAU  =  0,  but  TAU  ■  8160,  or  day  340.  The  time  of  year  is 
determined  by  a  variable  (SDEDY)  within  the  C-array.  For  experiments 
at  RAND,  day  340  corresponds  to  December  1. 

Pre-C0MP3  Arrays 

Rather  than  save  on  disk  all  the  variables  which  may  be  needed, 
the  values  of  some  state  variables  are  saved  before  they  are  adjusted 
by  subroutine  C0MP3,  thereby  allowing  the  variables  defined  in  COMP3 
to  be  obtained  by  a  postprocessing  pass.  The  pre-COMP3  variables  saved 


are : 


Table  3 


PRF-<'0MP3  ARRAYS 

T  Temperature  at  tropospheric  levels  1 
and  3  in  degrees  Kelvin 
Q3  Mixing  ratio  at  level  3 
GT  Ground  temperature  in  degrees  Kelvin 
GW  Ground  wetness  fraction 

State  Variables 

The  state  variables  after  the  complete  integration  step  are  stored 
on  the  disk  for  each  output  time.  The  arrays  are  dimensioned  46  x  72, 

representing  the  latitude/longitude  grid  points.  The  state  variables 
include : 


Table  4 

STATE  VARIABLES 

P  Surface  pressure  in  millibars,  minus  the  troposphere 
pressure  (200  mb) 

U  Horizontal  zonal  wind  at  levels  1  and  3  in  meters/sec 
V  Meridional  wind  at  levels  1  and  3  in  meters/sec 
T  Temperature  at  tropospheric  levels  1  and  3  in  degrees 
Kelvin 

Q3  Mixing  ratio  at  level  3 
GW  Ground  wetness  fraction 
GT  Ground  temperature  in  degrees  Kelvin 


For  further  explanation  of  these  symbols  and  arrays,  see  Ref.  1. 
Precipitation  Arrays 

Since  precipitation  rates  can  be  highly  variable  and  thus  the  true 
total  precipitation  be  lost  by  sampling  only  every  six  hours,  the  M/A 
program  accumulates  the  precipitation  data  every  half-hour.  This  in¬ 
formation  is  written  onto  dirk  in  three  arrays  for  each  six  hours  of 
simulation . 
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Table  5 

PRECIPITATION  ARRAYS 


1.  4CPPLS,  representing  PREC  *  P,  where  PREC  is 
level  3  moisture  change  due  tc  large-scale  pre¬ 
cipitation,  and  P  is  surface  pressure.  To  ob¬ 
tain  units  of  mm/day,  multiply  by  200/GRAV, 
where  GRAV  =  9.81  m/sec2. 

2.  ACPPCl,  representing  (PCI  +  Cl)  *  P,  where  PCI 
is  level  1  temperature  change  due  to  penetrating 
convection,  and  Cl  is  level  1  temperature  change 
due  to  mid-level  convective  latent  heating.  To 
obtain  units  of  mm/day,  multiply  by  200/ (GRAV  * 

CLH)  where  GRAV  =9.81  and  CLH  =  .24/580. 

3.  ACPPC3,  representing  (PC3  +  C3)  *  P,  the  equiv¬ 
alent  values  for  level  3.  To  obtain  units  of 
mm/day,  multiply  by  200/GRAV  *  CLH)  where  GRAV  = 

9.81  and  CLH  =  .24/580. 

ACCUMUiATED  ZONAL  MEANS  OUTPUT 

At  the  same  time  that  the  precipitation  data  set  is  written,  the 
M/A  program  writes  into  another  data  set  a  set  of  twenty-four  or  twenty- 
five  accumulated  zonal  arrays,  depending  on  the  version  of  the  M/A  pro¬ 
gram  used.  The  zonal  arrays  are  accumulated  every  half-hour  and  are 
written  onto  the  disk  every  six  hours.  At  each  time  point,  TAU ,  the 
time  in  hours,  is  written  in  one  record,  followed  by  the  twenty-four 
or  -five  records,  each  containing  a  46-word  zonal  mean  array. 

To  insure  that  the  zonal  arrays  are  stored  onto  disk,  the  M/A 

model  writes  TAU  and  the  zonal  arrays  followed  by  a  dummy  record  at 

each  time  point.  Then,  before  the  next  write,  a  backspace  over  the 

dummy  record  is  performed.  But,  at  the  end  of  the  data  set,  the  dummy 

record  remains.  This  record  can  be  recognized  by  the  fact  that  the 

—82 

first  value  is  much  smaller  than  TAU,  approximately  10  .  When  two 

of  these  data  sets  are  concatenated,  TAU  mus^  be  tested  so  that  the 
dummy  record  can  be  skipped. 

As  the  disk  is  written  every  six  hours,  and  the  computational  sub¬ 
routine  C0MP3  which  computes  tnese  values  is  executed  every  half-hour, 
each  array  on  disk  is  the  sum  of  the  arrays  for  12  time  points.  There¬ 
fore,  the  set  of  data  arrays  written  at  each  time  point  is  the  accumu¬ 
lation  from  the  previous  six  hours. 
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The  array  values  are  described  below  in  the  order  in  which  they 
appear  on  the  disk.  A  short  mnemonic  and  a  description  is  given  for 
each  variable.  If  a  formula  is  given,  the  names  used  in  th.  formula 
are  described  below  or  in  the  FORTRAN  dictionary  in  Ref.  1.  If  a  form¬ 
ula  is  not  given,  the  mnemonic  itself  is  a  variable  name  described  in 
that  reference. 

The  following  variable  names  are  used  in  the  formulas: 

L  =  latent  heat  of  condensation  =  580  cal/g 
CP  =  dry  air  specific  heat  at  constant  temperature  =  0.24 
cal/ g/ deg 

2 

g  =  gravitational  constant  =  9.81  m/sec 
p  =  surface  pressure  values,  which  change  with  time 

The  Zonal  Array  Values  in  disk  order  are  given  in  Table  6. 

Table  6 

ZONAL  ARRAY  VALUES  IN  DISK  ORDER 

1.  CL1 — the  deep  convective  cloudiness  indicator;  see  the 
Map  14  description  in  the  MAPGEN  User's  Guide. 

2.  CL2 — the  large-scale  convective  cloudiness  indicator; 
see  the  Map  14  description  in  the  MAPGEN  User's  Guide. 

3.  CL3 — the  low-level  convective  cloudiness  indicator; 
c.ee  the  Map  14  description  in  the  MAPGEN  User's  Guide. 

4.  CL — the  total  cloudiness  indicator;  see  the  Map  14 
description  in  the  MAPGEN  User's  Guide. 

5.  E4 — the  surface  evaporation  rate. 

6.  F4 — the  surface  sensible  heat  flux. 

7.  AS1 — the  insolation  absorbed  by  the  upper  layer. 

8.  AS3 — the  insolation  absorbed  by  the  lower  layer. 

9.  R0 — the  net  upward  long-wave  radiation  flux  at  trop- 
opause . 

10.  R2 — the  net  upward  long-wave  radiation  flux  at  level 

2. 

11.  P.4 — the  net  upward  long-wave  radiation  flux  at  the 
surface . 
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n.  S4— the  total  flux  of  short-wave  radiation  absorbed  by 
the  ground. 

13.  RETOT— the  total  solar  radiation  reflected  from  the 

Earth-atmosphere  system.  The  variable  RETOT  is  calcu¬ 
lated  in  the  following  way: 

RETOT  =  SS  *(-CIj  —  ALAC  *  (l  ~  ALS)  (1  -  CL)  *  ALAQ  *  (1  -  ALS) 

V  1  "  ALAC  *  ALS  1  -  ALAO  *  ALS 

+  FS4C  +/S4  *  ALS  TEMS  *  CL  \ 

\  1  -  ALS  1  -  TEMP  *  ALS/ 

where  TEMS  and  TEMP  are  one  of  four  sets  of  values  depending  on  the 
cloudiness  values  CL1,  CL2,  and  CL3,  as  follows: 

1.  If  CL1,  CL2,  and  CL3  s  0.0,  TEMS  =  TEJfP  =  0.0 

2.  If  CL3  :•  0.0  and  CL1,  CL2  s:  0.0,  then: 

TEMS  =  SA  *  (1  -  ALC3)  *  (1  -  .271  *  [ (EFV0  -  EFV3)/COSZ 
+  1.66  *  (EFVC3  +  EFV3)]’303) 

TEMP  =  ALC3 

3.  If  CL2  >  0  and  CL1  £  0,  then: 

TEMS  =  SA  *  (1  -  ALC2)  *  (1  -  .271  *  [ (EFV0  -  EFV2)/C0SZ 
+  1.66  *  EFVC2  *  COT  EFV3]’303) 

TEMP  =  ALC2 

where  COT  =  Min(1.0,  (RH3  -  0.9)  *  10.0) 

4.  If  CL1  >  0,  then: 

TEMS  =  SA  *  (1  -  ALC1)  *  (1  -  .271  *  [(EFV0  -  EFV1)/C0SZ 
+  1.66  *  (EFVC1  +  EFV3)]'303) 


TEMP  =  ALC1 . 
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14.  SCOSZ — the  total  solar  radiation  at  the  top  of  the  atmosphere. 

15.  STRA  -  SA  -  ASOT — the  short-wave  absorption  above  the  tropo¬ 

sphere. 

16.  PPTOT  -  [PREC  +(CP/L)  *  (Cl  +  C3  +  PCI  +  PC3) ]  (P/2g)  *  4800— 

the  total  precipitation. 

17.  HL1  -  (Cl  +  PCI)  *  48 — latent  heating  at  sigma  level  1/4. 

18.  HL3  =  (C 3  +  PC3  +  (CP/L)  *  (PREC)  *  48— 3  atent  heating  at 

sigma  level  3/4. 

19.  HSW1  -  AS1  *  (2g/P)  /  (CP/10)— the  short-wave  heating  at  sigma 

level  1/4. 

20.  HSW3  -  AS 3  *  (2g/P)  /  (CP/10)— the  short-wave  heating  at  sigma 

level  3/4. 

21.  HLW1  -  (R2  -  R0)  *  (2g/P)  /  (CP/10)— the  long-wave  heating  at 

sigma  level  1/4. 

22.  HLW3  -  (R4  -  R2',  *  (2g/P)  /  (CP/10)— the  long-wave  heating  at 

sigma  level  3/4. 

23.  HF4  =  F4  *  (2g/P)  /  (CP/10)— the  sensible  heat  from  the  Eartl  . 

24.  CICE  -  the  heat  rising  through  ice;  see  the  Map  23  description 

in  the  MAPGEN  User's  Guide. 

25.  HD1F — the  smoothed  heating  terms  as  computed  in  COMP 3 .  HI 

and  H3  are  the  temperature  changes  at  sigma  surfaces 
1/4  and  3/4  respectively.  These  values  are  smoothed 
before  being  added  to  the  temperature,  by  area  smoothing 
(HI  +  H3)/2  and  adding  (Hl-H3)/2  for  sigma  level  1/4  and 
subtracting  (Hl-H3)/2  for  sigma  level  3/4  temperature 
changes.  HDIF  is  the  smoothed  value  of  (Hi  +  H3)/2. 

Note:  This  zonal  array  was  not  put  onto  disk  in  the  version  of  the  M/A 
program  used  in  experiments  10-13.  The  error  in  calculating  total 
heating  with  this  term  omitted  is  not  large;  however,  the  array  has 
been  added  for  the  sake  of  accuracy. 


The  values  saved  on  disk  differ  by  a  constant  from  the  six-hour 
average  quantities  listed  above.  Each  of  these  arrays  has  to  be  multi¬ 
plied  by  the  constant  given  in  Table  7  in  order  to  obtain  desired  units. 
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as  shown.  The  number  12  appears  in  the  constants  because  C0MP3  is  ex¬ 
ecuted  12  times  in  each  six-hour  period,  and  72  is  the  number  of  longi¬ 
tude  points  accumulated  at  each  latitude.  The  radiation  unit  of  lang— 
leys  is  abbreviated  as  ly. 
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Table  7 

CONVERSION  FACTORS  FOR  ZONAL  ARRAYS  ON  DISK 


Variable  x  Const.  =  Units 

Variable 

>  x  Const.  = 

Units 

1 

1 

CL1 

12*72 

fraction 

RETOT 

12*72 

ly/day 

1 

1 

CL2 

12*72 

fraction 

SCOSZ 

12*72 

ly/day 

1 

1 

CL3 

12*72 

fraction 

STRA 

12*72 

ly / day 

1 

100 

CL 

12*72 

fraction 

PPTOT 

36  *g 

mm/day 

48 

E4 

1000 

mra/day 

HL1 

12*72 

deg  C/day 

1 

48 

F4 

12*72 

ly/day 

HL3 

12*72 

deg  C/day 

1 

2*8 

ASl 

12*72 

ly/day 

HSW1 

12*72* 10*CP 

deg  C/day 

1 

2*8 

AS3 

12*72 

ly/day 

HSW3 

12  •  72  *10  *CP 

deg  C/day 

1 

2*8 

RO 

12*72 

ly/day 

HLW1 

12*72* 10*CP 

deg  C/day 

1 

2*8 

R2 

12*72 

ly/day 

HLW3 

12  *72  *10*CP 

deg  C/day 

1 

2*g 

R4 

12*72 

ly/day 

HF4 

12 • 72  *10  *CP 

deg  C/day 

1 

1 

S4 

12*72 

ly/day 

CICE 

12*72 

ly/day 

48 

HDIF 

12*72 

deg  C/day 
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III.  USER'S  GUTDE  FOR  MAPGEN 


The  MAPGEN  program  reads  the  disk  written  during  the  M/A  simula¬ 
tion  program,  combines  variables  stored  on  :he  disk  into  meteorologi¬ 
cally  meaningful  parameters,  and  prepares  a  tape  for  the  Map  Display 
program  which  draws  contour  maps  and  graphs  of  the  zonal  averages.  The 
standard  output  (see  description  below)  consists  of  maps  and  zonal  av- 
erage  graphs  of  28  variables.  Many  of  the  variables  can  be  mapped  at 
more  than  one  level.  While  the  28  variables  are  commonly  used,  new 
maps  may  be  defined  and  included  in  MAPGEN,  or  some  of  the  28  may  be 
excluded . 

In  order  to  recover  the  meteorological  variables  calculated  during 
the  integration  of  the  M/A  model  within  subroutine  COMP 3 ,  MAPGEN  uses 
the  pre-C0MP3  arrays  on  disk  plus  the  subroutine  COMP 3  to  recalculate 
these  values.  The  version  of  COMP 3  used  for  MAPGEN  is  different  from 
the  M/A  COMP 3 ,  but  the  actual  calculations  are  the  same  for  both  pro¬ 
grams.  Many  experiments  involve  a  change  to  COMP 3  in  the  M/A  program, 
and  these  changes  must  also  be  made  to  the  MAPGEN  version  of  COMP 3  in 
order  to  obtain  correct  values.  For  this  reason,  as  a  general  rule, 
COMP 3  should  always  be  compiled  when  MAPGEN  is  run.  This  listing  can 
then  be  referred  to  if  questions  arise  about  the  values  calculated  in 
MAPGEN. 

MAPGEN  can  sample  data  from  the  history  disk  at  specified  single 
time  points  or  collect  and  time-average  data  over  specified  time  per¬ 
iods.  In  calculating  the  time  averages,  the  starting  and  ending  data 
points  are  given  half  the  weight  of  each  interior  time  point.  If  these 
times  were,  for  example,  370.0  to  400.0 — a  thirty-day  average — the 
values  at  each  end  point  would  be  multiplied  by  one-half  and  then  added 
to  the  values  for  each  six  hours  from  370.25  to  399.75  days.  At  the 
end  of  the  accumulations,  the  numbers  are  averaged  by  dividing  the 
total  by  120 — not  by  121,  as  would  be  the  case  if  the  beginning  and 
ending  days  were  counted  at  full  value.  Although  a  thirty-day  average 
is  most  commonly  requested,  the  averaging  period  is  arbitrary  and, 
like  the  other  options,  is  controlled  by  tne  input  cards. 
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INPUT  FORMATS 

The  first  card  of  the  MAPGEN  run  Is  the  header  card  (read  only 
once  for  each  run). 


HEADER  CARD  FORMAT 


Co] jmn 

Format 

Meaning 

1 

11 

Describes  new  disk  format;  leave 

blank 

2-77 

19A4 

^  title  printed  at  the  beginning 
of  MAPGEN  output.  In  hap  Display, 
this  title  is  printed  on  each  page 
and  frame  of  output. 

Following  the  header  come  the  sets  of  map  cards.  The  first  card 
of  each  set  inputs,  among  other  things,  the  start  and  end  times  for 
the  time  averaging.  In  MAPGEN,  a  maximum  of  40  maps  may  be  called  for 
with  one  set  of  data;  however,  in  each  run,  as  many  sets  may  be  input 
as  needed.  Each  set,  if  it  contains  less  than  40  map  cards,  is  term¬ 
inated  by  a  blank  card,  and  the  last  set  is  terminated  by  two  blank 

cards.  (If  exactly  40  map  cards  are  input,  the  sets  should  not  be 
separated  by  a  blank  card.) 

If  a  run  consists  of  more  than  one  set  of  data,  the  user  should 
keep  in  mind  that  the  computational  subroutine,  C0MP3,  uses  most  of 
the  CPU  time.  Thi*  is  true  even  though  all  of  the  variables  computed 
by  COMP 3  at  a  time  point  are  found  in  one  pass  through  the  subroutine. 
The  13  variables  that  use  COMP 3  should  therefore  be  collected  into  one 
set  of  data,  and  those  that  do  not,  into  another.  The  map  numbers  of 
the  15  variables  that  do  not  use  COMP 3  are  1-8,  10-13,  and  15-17.  The 

format  of  the  first  card  of  each  set  is  given  in  Table  8,  and  the  for¬ 
mat  for  the  map  cards,  in  Table  9. 

The  data  cards  for  MAPGEN  can  also  be  used  for  the  Map  Display 
program,  which  minimizes  input  preparation  for  the  two  programs.  How¬ 
ever,  for  Map  Display,  a  blank  card  must  always  separate  the 
input  cards. 


sets  of 
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Table  8 

FORMAT  OF  FIRST  CARD  OF  EACH  SET  FOR  MAPGEN  INPUT 


Column 

Format 

Meaning 

1-10 

F10.0 

Time  to  start  in  days,  e.g.,  370.0 

11-20 

F10.0 

Time  to  end  in  days,  e.g.,  400.0 

Note:  The  data  for  a  single  time 
point  can  be  obtained  by  setting 
the  start  and  end  times  to  be  the 

same . 

21-30 

F10.0 

Projection  type: 

0.0  =  rectangular 

1.0  =  north  stereographic 
-1.0  =  south  stereo 

2.0  =  north  equidistant 
-2.0  =  south  equidistant. 

31-40 

F10.0 

If  nonzero,  cards  of  the  zonal  mean 
values  are  punched  in  flAPGEN.  The 
value  in  these  columns  is  truncated 
to  an  integer  value  and  punched  on 
the  cards  as  the  experiment  number. 
See  Punched  Output  described  on 
page  23. 

I 
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T ab  It*  9 

KORMAT  OK  MAI*  CARDS  FOR  MAI’CLN  AND  MAI*  DISPLAY  INPUT 


Column 

Format 

Meaning 

1-2 

12 

The  map  number,  between  1  and  28,  right 

adjusted . 

1 

F5.0 

SURF,  the  sigma  surface,  e.g.,  0.25,  0.75, 

1.0,  or  for  some  maps  the  pressure,  e.g., 

400.0  mb,  800.0  mb.  For  convenience,  this 

value  will  be  called  SURF  in  later  discus¬ 
sions.  "'he  sigma  surface  is  1.0  for  the 

surface  and  0.0  for  the  top  of  the  tropo¬ 
sphere. 

F5.0 

The  isoline  interval.  Isoline  values  are 

calculated  using  positive  and  negative 
multiples  of  this  value,  starting  from  the 

dashed-line  isolin'-  value. 

13-17 

F5.0 

Value  of  the  isoline  tn  be  dashed. 

18-22 

F5.0 

Not  used. 

23-27 

F5.0 

Number  of  copies  of  the  array  to  be  printed 
at  this  time,  normally  1.0.  If  0.r,  the 

values  are  not  printed. 

28-32 

F5.0 

Number  of  specific  isoline  values  to  be  read 
(£15).  If  these  columns  are  zero  or  blank, 
the  isoline  intervals  described  by  8-12  and 
13-17  are  used.  Otherwise,  one  additional 
card  is  read  for  the  map  description.  It 

contains  the  map  number  in  columns  1-2  as 

usual,  followed  by  as  many  as  15  specific 

values  of  isolines  to  be  contoured  (each  tn 

F5.0  format).  This  extra  card  d~>es  not  count 
toward  the  limit  of  40  map  cards  per  sc it. 

Note  that  the  dash  line  value  has  to  be  in¬ 
put  among  the  values  if  it  is  to  be  drawn. 

33-37 

F5.0 

If  not  zero,  then  zonal  means  are  not  to 

be  graphed  in  the  Map  Display  program. 

38-42 

F5.0 

If  not  zero,  the  contour  map  is  not  graphed 

in  the  Map  Display  program. 
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Only  the  first  seven  columns  of  the  map  cards  are  used  by  MAPGEN . 
The  map  number  determines  the  variable  to  be  used,  and  the  value  of 
SURF  specifies  the  level  at  which  the  variable  is  desired. 

STANDARD  OUTPUT  MAP  DESCRIPTIONS 

A  standard  set  of  variables  is  plotted  for  each  simulation.  This 
standard  output  package  consists  of  two  sets  of  data,  the  first  not  us¬ 
ing  C0MP3  and  the  second  using  the  subroutine .  The  card  input  for  this 
package  is  listed  in  Appendix  A.  Some  of  the  28  different  variables  are 
available  at  more  than  one  level.  The  val ue  of  SURF  (meaning  surface 
or  pressure)  usually  denotes  the  sigma  surface  level  but  it  is  some¬ 
times  a  pressure  level  or,  in  some  cases,  simply  a  flag  to  distinguish 
different  maps. 

Table  10  gives  the  map  numbers  and  titles,  showing  the  units  of 
the  variables.  The  values  of  SURF  are  also  given  for  the  standard  in¬ 
put  package — SURF  usually  indicates  the  tropospheric  level,  expressed 
either  as  a  sigma  level  or  as  a  pressure  in  millibars,  as  described 
above.  Some  of  the  variables  are  identified  by  different  map  numbers 
in  Chap.  IV  of  Ref.  1,  the  Mintz-Arakawa  documentation.  The  variables 
have  since  been  renumbered.  See  Appendix  C  for  the  correspondence  of 
map  numbers.  Variables  not  discussed  in  Ref.  1  are  described  in  Table 
11,  below. 

MAPGEN  OUTPUT 

The  printed  output  of  MAPGEN  is  very  limited.  The  title  card  is 
printed,  then  the  input  cards  of  each  data  set  are  printed  out,  with 
the  end  of  each  set  followed  by  a  row  of  zeros.  After  the  data  set 
has  been  printed,  the  M/A  disk  data  sets  are  read  in.  For  each  time 
record  read,  the  time  point  in  hours  and  days  is  written  out.  A  new 
page  is  started  when  the  first  time  point  to  be  included  is  reached. 
When  all  of  the  values  have  been  averaged,  the  information  is  written 
onto  tape  and  a  description  of  each  map  is  written  out,  consisting  ot 
the  first  and  last  times  of  the  accumulation,  the  number  of  accumula¬ 
tions,  the  map  number,  the  SURF  value,  and  the  map  title. 
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Table  10 

MAP  TITLES  FOR  MAPGEN  OUTPUT 


Map 

No. 

SURF  fo ;  Standard 

Output 

Uses 
COMP  3 

Title 

1 

1.0 

No 

Topography  (ocean  temperature,  deg  C) 

3.0 

No 

Topography  (surface  elevation,  hectofeet) 

2 

1.0 

No 

Sea  Level  Pressure  Smoothed  (mb  -  1000) 

3 

0.25,  0.75,  1.0 

No 

Pressure  at  Sigma  Surface  (mb) 

4 

0.25,  0.75,  400,  800 

No 

Geopotential  Height  (hectometers) 

5 

0.25,  0.75,  400  800 

No 

East-West  (u)  Wind  Component  (m/sec) 

6 

0.25,  0.75,  400,  800 

No 

North-South  (v)  Wind  Component  (m/sec) 

7 

0.5 

No 

Sigma  Vertical  Velocity  (mb/hr) 

8 

1.0 

No 

Ground  Temperature  (deg  C) 

9 

1.0 

Yes 

Surface  Air  Temperature  (deg  C) 

10 

0.25,  0.75,  400,  800 

No 

Temperature  (deg  C) 

11 

1.0 

No 

Ground  Wetness  (scaled  zero  to  ten) 

12 

0.75 

No 

Relative  Humidity  (percent) 

13 

1.0 

No 

Precipitable  Water  (cm) 

0.25 

Yes 

CL1  -  Deep  Convective  Cloudiness 

14 

0.75 

Yes 

CL2  -  Large-Scale  Cloudiness 

1.0 

Yes 

CL3  -  Low-Level  Convective  Cloudiness 

2.0 

Yes 

CL  -  Total  Cloudiness 

15 

1.0 

Nc 

Large-Scale  Precipitation  (mm/day) 

16 

1.0 

No 

Convective  Precipitation  (mm/day) 

17 

1.0 

No 

Total  Precipitation  (mm/day) 

lb 

1.0 

Yes 

Evaporation  (mm/day) 

19 

1.0 

Yes 

Surface  Heating  by  Evaporation  (100  ly/day) 

20 

1.0 

Yes 

Sensible  Heat  Flux  (10  ly/day) 

21 

1.0 

Yes 

Surface  Insolation  Absorption  (100  ly/day) 

22 

1.0 

Yes 

Surface  Long-Wave  Heating  (100  ly/day) 

23 

1.0 

Yes 

Surface  Heat  B  ilance  (100  ly/day) 

24 

0.25 

Yes 

Latent  Heating  in  Layers  (sigma  ■  0  to  . 5)(deg'day) 

0.75 

Yes 

Latent  Heating  in  Layers  (sigma  -  .5-1)  (deg/day) 

25 

0.25,  0.75 

Yes 

Absorption  of  Insolation  in  Layers  (deg  C/day) 

26 

0.25,  0.75 

Yes 

Long-Wave  Heating  in  Layers  (deg  C/day) 

27 

0.25,  0.75 

Yes 

Total  Heating  (deg  C/day) 

28 

U _ 

1.0 

Yes 

Net  Albedo 
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Table  11 

VARIABLES  NOT  COVERED  IN  GATES ,  ET  AL .  ,  BY  MAP  NUMBER 


Map 

Variable 

Formula  or  Calculation  Method 

3 

Pressure 

PTROP  +  SURF  •  3 

where  PTROP  *  tropopause  pressure  -  200mh ,  SURF  -  sigma  surface,  P  -  surface 
pressure  parsmeter 

7 

Sigma  Vertical  Velocity 
Meridional  Mass  Flux 

Defined  in  Ref.  1,  p.  170.  From  these  values,  the  zonal  mean  values  of  the  merid¬ 
ional  mass  flux  are  calculated;  these  sre  printed  and  graphed  in  the  Map  Display 
program.  Calculated  at  sigma  surface  of  1/2.  SURF  is  not  used.  See  p.  28  for 
tormula . 

U 

C loudtness 

CL1  Deep  convection 

CL2  Large-scrle 

CL3  Low- level 

CL  Total 

CL1 

CL2 

CL  1 

CL 

Collected  when  SURF  0.5 

Collected  when  SURF  between  0.5  and  1.0 

Collected  when  SURF  over  1.0 

Calculated  in  COMP3  by: 

Deep  Convective  Cloudiness  (CM) 

if  both  Cl  and  PCI  70.0,  CL1  ■  0.0,  i.e.,  if  there  is  no  positive  temperature  change 
in  level  1  due  to  convective  latent  heating  (Cl)  or  to  penetrating  convection  (PCI); 
otherwl se  , 

CL1  -  Max  [0.0,  Min (RH3  •  1.25)-. 25,  1.0)] 
where  RH  3  is  the  level  3  relative  humidity. 

Large-Scale  Cloudiness  (CL2) 

CL2  Is  either  0.0  or  1.0.  CL2  -  1.0  If  RHT  0.9,  oi heru.se  CL2  -  0.0. 

Low-Level  Convective  Cloudiness  T CL 3) 

CLT  Is  0.0,  except  If  the  stability  parameter  for  low-1  vel  convection  (EX)  Is  greater 
than  zero  and  the  level  1  temperature  change  due  to  p  i  trstlng  convection  (PCI)  Is 
equal  to  0.0,  In  which  case: 

CLT  -  Max(0.0,  Min (1 . 25  •  RHT)-.:  ,  1.0)1 

Total  Cloudiness  (CL) 

CL  -  Max  (Cl.l,  CL2 ,  Cl.)) 

For  these  maps,  Sl'RF  Is  not  printed  on  the  map  frames. 

15 

Large-scale  Precipitation 

Computed  from  the  accumulated  half-hour  values  of  PRFC  •  P  (precipitation  at  each  grid 
point)  that  are  saved  at  each  output  time.  It  is  given  by  the  expression: 

48  102  (  1  V"  (PREC  ■  P)  \ 

2* '  v2  h  ) 

where  "1  g/cm,  i  represents  the  ith  half-hour  value,  PRFC  is  .he  level  )  moisture 

w 

change  due  to  large-scale  condensation,  and  P  is  the  surface  pressure  parameter.  SURF 
is  not  used. 

16 

Convective  Precipitation 

Produced  from  the  accumulated  half-hour  values  of  (PCI  +  Cl)  •  P  and  (PCI  +  Cl)  •  " 

that  are  saved  at  each  output  time.  It  is  given  bv  the  expression: 


48  . 

2g  °v 

SURF  is  not  used. 


' 

12 

12 

' 

i 

£  l (PCI  +  Cl) 

P]  +  5Z  [  (PCT  +  C3)  •  P] 
1-1 

1-1 

Formula  or  Calculation  Method 


Surface  Heating  by 
Evaporation 


Surface  Long-vave 
Heat lng 


Surface  Heat  Balance 


Sigma  0. 5-1.0 


Computed  from  the  expresalon: 


-L  •  E4  •  Rh4.0 


Latent  Heating  In  Layers 


Sigma  0.5 


Where  L  latent  heat  (580.0)  and  E4  -  surface  evaporation  rate.  SURF  la  not  used. 

The  negative  of  this  variable  Is  described  on  p.  212  of  Ref.  1.  The  expression  Is: 

-R4/100 

where  R4  la  the  net  upward  long-wave  radiation  at  the  Earth's  surface.  SURF  la  not 


The  values  for  this  -nap  equal 

Map  21  +  Map  22  -  (Map  20  •  10  *)  +  Map  19  +  C1CF.  •  10_^ 
where  C1CE  Is  0  for  land  or  Ice-free  ocean.  When  this  grid  point  Is  sea  Ice, 

CICE  -  TEM  •  (TICE  -  TC) 

thtJknir  r  ***  I1"""*’  conductlvlty  °f  Ice  (-432'  over  the  effective  Ice 

tuJe  SURF  is  n^t  used'  Mlt1"8  P°lnt  °f  9"d  TC  19  th*  *r°“"d 


The  map  for  the  sigma  level  0  to  one-half  Is  described  by: 

(Cl  +  PCI)  •  48 

where  Cl  and  PCI  are  level  1  temperature  change  due  ty  midjevel  convective  latent 
heating  and  to  penetrating  convection,  respectively.  The  number  48  puts  the  values 
in  the  units  degs  C/dsy. 

The  map  for  the  algma  level  one-half  to  1.0  Is  described  by: 

(C3  +  PC3  +  CL1I  •  PRF.C)  •  48 

|*ere  C3  and  PC3  are  level  3  temperature  change  due  to  mid-level  convective  latent 
heating  and  to  penetrating  convection,  respectively,  CLH  Is  the  latent  heat  to  spe¬ 
cific  heat  ratio  (580/. 24),  and  PREC  la  level  3  moisture  change  due  to  large-scale 


Computed  from  the  ratio  of  two  time-averaged  variables.  At  each  grid  point,  the  time- 
averaged  values  of  the  following  variables  are  collected: 

SCOSZ  and  RET0T 

where  RET0T  Is  the  total  reflected  solar  radiation,  and  SCOR2  Is  the  total  solar  rad¬ 
iation  at  the  top  of  the  atmosphere.  (See  Sec.  II  for  the  formulation  of  RETOT.)  The 
variable  values  are  then  formed  using  RET0T/SC0R2.  SURF  la  not  u'-ed 
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Tape  Output 

The  important  output  of  MAPGEN  is  the  tape  that  contains  the  map 
data.  For  each  map,  a  tape  record  is  written  describing  the  map  data. 
These  records  are  used  as  input  to  the  Map  Display  Program,  where  the 
actual  contour  and  zonal  profile  maps  are  drawn.  If  two  experiments 
are  to  be  compared,  the  data  arrays  can  be  input  to  the  differencing 
program,  which  subtracts  them,  and  Map  Display  can  contour  the  differ¬ 
ences  . 

The  following  information  is  put  into  each  record  by  MAPGEN. 


1. 

2. 


3. 

4. 

5. 

6. 
7. 


8. 

9. 

10. 

11. 

12. 


The  start  and  end  times  in  hours  of  the  accumulation  (2  words) 
The  (floating-point)  weighting  factor  used  to  average  the 
accumulations.  If  only  one  time  point  is  used,  this  number 
is  set  to  1.0. 

A  4-character  identification,  input  in  the  M/A  run. 

The  map  number  (integer). 

The  52-character  13-word  title  of  the  map. 

The  sigma  surface  or  pressure,  SURF. 

Two  logical  flags,  one  word  each,  stating  whether  the  array 
data  are  on  the  n  grid  (the  horizontal  finite-difference 
grid)  or  the  u,v  wind  velocity  grid.  See  Ref.  1,  p.  70. 

A  word  no  longer  used. 

The  46  x  72  array  values. 

The  46-word  zonal  mean  of  the  variable. 

The  46-word  zonal  mean  of  the  variable  over  the  oceans. 

The  global  mean. 


A  dummy  record  of  the  same  length  is  used  to  flag  the  end  of  the  tape. 
The  start  and  end  times  for  this  record  are  set  to  values  over  100,000. 


Punched  Output 

If  the  punched  output  for  the  zonal  mean  values  is  requested  by 
the  first  card  of  the  data  set,  then  for  each  map  in  the  data  set, 
seven  punched  cards  giving  the  zonal  means  and  the  global  mean  will 
be  generated.  The  first  three  columns  of  the  cards  give  the  card 
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number,  0  through  6.  The  next  six  columns  give  the  input  (integer)  ex¬ 
periment  number  and  the  map  number,  followed  by  the  SURF  value  in  col¬ 
umns  10-16.  The  first  card  gives  the  start  and  end  times  in  columns 
17-24  and  25-32.  The  next  five  cards  have  eight  zonal  mean  values  per 
card  in  columns  17  through  80,  8  columns  per  value,  with  the  values 
given  to  thousandths.  The  last  card  has  the  last  six  zonal  mean  values 
plus  the  global  mean  in  columns  65-72. 

JOB  CONTROL  COMMENTS 

The  program  is  run  at  UCLA  on  the  IBM  360/91.  The  program  uses 
1250K  of  core,  and  is  compiled  on  the  H-compiler.  The  map  data,  along 
with  other  information,  is  written  onto  FORTRAN  unit  9,  in  this  case 
a  tape,  in  one  large  record  for  each  map. 

The  disk  data  sets  written  by  the  M/A  program  are  input  on  FORTRAN 
unit  8.  Two  disk  data  sets  may  have  to  be  concatenated  as  one  disk 
holds  only  about  23  days  simulation  time. 
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IV.  USER'S  GUIDE  FOR  MAP  LISPLAY 


Map  Display  outputs  the  data  generated  by  MAPGEN.  The  data  can  be 
displayed  in  three  ways:  The  array  and  affiliated  statistics  can  be 
printed,  a  contour  map  can  be  drawn,  or  the  zonal  mean  values  can  be 
graphed.  The  type  or  types  of  output  are  controlled  by  card  input. 

The  choice  of  the  isolines  to  be  drawn  in  the  contour  map  is  very 
important.  These  values  may  be  a  personal  choice,  or  may  depend  on 
the  data;  one  experiment  may  need  different  isolines  from  another. 

This  is  one  reason  why  Map  Display  and  MAPGEN  are  two  programs;  the 
array  data  may  be  used  to  generate  many  different  contour  maps  with 
different  isolines.  However,  the  isolines  used  in  the  standard  output 
package  are  very  good  for  most  runs.  (See  Appendix  A.) 

INPUT  DESCRIPTION 

There  are  three  types  of  input  for  Map  Display:  the  map  descrip¬ 
tions  similar  to  the  MAPGEN  input,  the  land  outlines  for  the  contour 
maps,  and  the  array  data  from  MAPGEN.  The  first  two  types  are  input 
with  cards,  and  the  MAPGEN  output  is  on  tape  or  disk. 

To  obtain  an  outline  of  the  continents  used  in  the  M/A  program, 
a  set  of  input  cards  must  be  prepared  for  the  Map  Display  program. 

For  each  M/A  grid  point,  a  \ alue  is  input:  zero  if  the  point  is  ocean, 
and  nonzero  if  the  point  is  land.  Note  that  sea  ice  is  considered 
ocean  (0)  and  land  ice  is  considered  land  (1) . 

The  land  mass  data  are  input  using  144  cards.  For  eacli  of  72 
meridians,  two  cards  are  used  to  describe  the  46  latitude  values  along 
the  meridian.  The  first  card  has  25  values  in  the  first  75  columns 
in  F3.0  format  (0.0  or  1.0),  starting  with  the  South  Pole,  and  the 
second  card  has  the  remaining  21  values  in  the  first  63  columns  in 
F3.0  format,  ending  at  the  North  Pole.  The  72  pairs  of  cards  start 
at  180°  W  longitude  and  proceed  eastward  to  175°  E. 

The  input  data  sets  described  in  the  MAPGEN  section  (pp.  16-19) 
are  also  used  for  Map  Display.  The  first  card  of  the  data  set  and  the 
map  descriptions  are  the  same  as  those  described  in  the  MAPGEN  User's 


-26- 


Guide.  However,  the  input  is  more  flexible  for  Map  Display.  For  one 
thing,  one  data  set  can  consist  of  any  number  of  map  descriptions, 
while  MAPGEN  is  limited  to  40.  Incidentally,  the  order  of  the  map 
descriptions  does  not  have  to  be  the  same  as  on  tape;  they  can  be  put 
into  any  order  desired.  To  match  the  data  on  tape,  only  the  start  and 
end  times,  the  map  number,  and  the  SURF  value  need  to  be  the  same;  all 
the  other  values  in  the  map  description  can  be  different  in  the  Map 
Display  run.  Thus,  in  one  Map  Display  run,  several  contour  maps  using 
different  isolines  can  be  drawn  using  the  same  array  data. 

PRINTED  OUTPUTS 

Appendix  B  presents  an  example  of  the  printout  from  the  Map  Dis¬ 
play  program.  Each  map  printout  consists  of  five  pages  of  printed  out¬ 
put,  with  a  standard  heading  on  each  page. 

The  output  heading  consists  of  three  lines  of  data.  The  first 
line  gives  the  map  number  and  the  run  title  that  were  input  in  columns 
2  through  77  of  the  first  card  of  the  input .  The  second  line  gives 
the  map  title,  the  SURF  value,  the  contour  interval,  and  rhe  global 
mean.  The  SURF  value  is  identified  as  S/P,  which  means  sigma  surface 
or  pressure  level,  and  the  contour  interval  is  called  CONINT.  The 
third  line  consists  of  the  four-character  identification  input  in  the 
M/A  run,  time  points  sampled  from  the  M/A  disk  by  MAPGEN,  and  the 
starting  and  ending  days. 

The  global  mean  is  calculated  by  area  weighting  the  zonal  means. 
There  are  two  sets  of  weighting  functions,  one  for  the  standard  it  lat¬ 
itude/longitude  grid  and  the  other  for  the  wind  grid  (which  is  called 
the  u,v  grid  in  the  M/A  model).  The  weighting  function  (wf)  for  the 
ir-grid  is  defined  by: 

(  1  -  Sin  88°  for  L  ■  -90°  and  +90° 

wf  (L)  -  . 

(  sin  (L  +  2°)  -  sin  (L  -  2°)  for  L  =  -86°  to  +86° 
where  L  proceeds  from  -90°  to  +90°  by  increments  of  4°. 
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The  weighting  function  for  the  u,v  grid  is  given  by: 

wf  (L)  =  sin  (L  +  2°)  -  sin  (L  -  2°)  for  L  =  -88°  to  +88° 

where  L  goes  from  -88  to  +88  by  4°  increments. 

The  global  average  of  a  zonal  mean,  X(L) ,  is  then  given  by: 

\  l  wf  (L)  •  X(L) 

The  latitudes  are  printed  on  the  left  side  of  each  page.  The  last 
row  on  the  page,  labeled  AVG,  represents  the  weighted  meridional  aver¬ 
ages  of  the  longitude.  The  longitude  labels  are  directly  under  the 
page  h  ading.  The  last  page  of  output  includes  the  zonal  means  and  the 
zonal  means  over  the  oceans. 

Whenever  a  local  high  or  low  value  is  found,  an  H  or  L  is  printed, 
respectively.  In  the  case  of  the  wind  maps  5  and  6,  these  letters  are 
replaced  by  W  or  E  for  map  5,  and  S  or  N  for  map  6. 

The  zonal  means  graph  has  the  global  mean  printed  at  the  right 
side  of  the  graph.  Below  this,  the  zonal  mean  values  are  printed  next 
to  their  respective  latitudes.  The  line  in  the  center  of  the  graph 
area  represents  the  zero  value.  The  scale  for  the  plot  is  determined 
by  the  zonal  means  themselves;  the  largest  absolute  value  among  the 
means  is  chosen  as  the  scale  factor.  In  the  example  in  Fig.  2,  the 
largest  absolute  value  is  -1.683,  and  thus  the  extreme  left  of  the 
graph  area  is  chosen  to  be  -1.683  and  the  extreme  right  of  the  graph 
area  is  +1.683.  The  hash  marks  on  the  zero  value  line  represent  each 
of  the  46  latitudes. 

For  map  7,  which  describes  the  sigma  vertical  velocity,  an  addi¬ 
tional  output  is  given.  Zonal  means  of  the  horizontal  mass  flux  of 
the  simulation  are  printed  with  the  other  means.  If  a  zonal  mean  graph 
for  map  7  is  to  be  plotted,  the  mass  flux  profile  will  be  drawn,  in  ad¬ 
dition  to  the  graph  of  vertical  velocity.  A  different  map  title  is 
used  for  the  horizontal  mass  flux  profile.  If  ZM(J)  is  the  original 
zonal  mean  then  ZMF(J),  the  zonal  mean  of  the  horizontal  mass  flux,  is 
defined  by: 
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ZMF(J) 


^0.0  if  J  =  1 

\  ZMF(J-l)  -  2.52  *  [ ZM ( J )  *  COS  (4°  *  (J-l)  -  86° 

(  +  ZM(J-l)  *  COS  (4°  *  ( J-2)  -  <16°)]  for  J  =  2 ,  46 


The  unit  of  measurement  is  10 


13 


gm/ sec . 


PLOT  OUTPUT 

Examples  of  the  two  plots  drawn  by  Map  Display  are  given  as  1  igs. 

1  and  2.  Figure  1  is  the  contour  map  and  Fig.  2  is  the  corresponding 
plot  of  zonal  means. 

The  heading  is  similar  to  the  printout  heading.  The  first  line 
consists  of  the  map  number  and  the  run  title.  The  second  line  starts 
with  the  map  title.  The  SURF  value  is  normally  printed,  either  with 
"SIG  ="  if  SURF  is  a  sigma  surface,  or  with  "MB"  following  the  value  if 
SURF  is  a  pressure  level.  For  maps  1  and  14,  SURF  is  not  printed  at 
all,  as  SURF  is  strictly  a  flag  in  these  maps.  Following  this,  on  the 
same  line,  is  the  message  ISOLINES  AT,  followed  by  the  smallest  isoline 
value  and  the  contour  interval.  The  third  line  starts  out  defining 
the  time  intervals,  giving  the  start  and  end  time  in  days.  Next  a  four- 
character  identification  is  printed,  which  was  input  during  the  M/A 
simulation  run.  Next  the  dashed-line  isoline  value  is  printed,  pre¬ 
ceded  by  the  message  DASHED  LINE  IS. 

For  maps  which  use  the  u,v  grid,  a  warning  message,  GRID  OFFSET, 
is  written  at  the  end  of  this  line.  That  is  because  the  contours  are 
drawn  as  if  the  values  are  on  the  ir-grid,  when  actually  the  contour 
lines  are  displaced  one-half  grid  interval  to  the  left  and  upward. 

The  contour  map  has  the  latitude  values  printed  every  20  degrees 
on  the  right  side  of  the  map.  At  each  grid  point  in  the  oceans,  a  dot 
is  plotted  so  that  the  latitude  and  longitude  are  more  easily  read. 

The  land  outline  is  drawn  between  the  land  and  water  grid  points.  Ice 
is  not  indicated  on  the  map. 

JOB  CONTROL  COMMENTS 

To  draw  the  maps,  Map  Display  uses  a  Rand-developed  software  sys¬ 
tem  called  the  Integrated  Graphics  Systems  (ICS) ,  described  in  Ref.  2. 
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MAP  17  EXP  DIFF.  EXP  16  MINUS  EXP  15  MAPS  13.  15.  16.  17 

TOTAL  PRECIPITATION  (MM  /  DAY)  SIG=  1. 000 

INTERVAL^  400.00  TO  430.00  LRL 


GLOBAL  r*AN>  -0. 190 


LAT 

JONAL  rCAN 

90. 

0.074 

66. 

-0.133 

86. 

-0.137 

78. 

0.137 

74. 

0.196 

70. 

-0.106 

66. 

-0.667 

66. 

0.054 

50. 

0.164 

54. 

0.663 

50. 

0.663 

46. 

-0.174 

46. 

-0.853 

38. 

-1.163 

34. 

-1.094 

30. 

-0.439 

66. 

0.169 

66. 

0.767 

18. 

1.571 

14. 

1.506 

10. 

0.104 

6. 

-0.701 

6. 

-1.683 

-6. 

-1.330 

-6. 

-1.063 

-10. 

-0.681 

-14. 

-0.695 

-18. 

-0.563 

-66. 

-0.019 

-66. 

0.359 

-30. 

0.017 

-34. 

0.103 

-38. 

-0.178 

-46. 

-0.654 

-46. 

-0.199 

-50. 

0.016 

-54. 

-0.606 

-58. 

-0.030 

-6 6. 

-0.134 

-66. 

-0.015 

-70. 

0.068 

-74. 

0.007 

-7E. 

0.016 

86. 

0.030 

-86. 

0.576 

-90. 

0.170 

Fig.  2 — Sample  plot  of  zonal  means  output  by  the  Map  Display 
The  accompanying  contour  map  is  shown  as  Fig.  1. 


program . 
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1GS  writes  a  7-trr  '*  ,  800  bpi  tape,  which  can  be  used  by  an  S-C  4060 
graphical  output  device  to  draw  the  pictures  on  microfilm  and  print 
them  on  11  x  14  paper. 

The  execution  step  uses  three  FORTRAN  units.  In  Rand's  operating 
system,  the  S-C  4060  output  uses  the  file  name  SC4060ZZ.  FORTRAN  unit 
8  is  the  output  MAPGEN-type  tape,  and  FORTRAN  unit  9  is  the  input  unit 
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V.  USER’S  GUIDE  FOR  THE  ZONAL  PROFILE  PROGRAM 


The  Zonal  Profile  program  graphs  the  time  averages  of  the  zonal 
arrays  written  by  the  M/A  program.  The  graphs  are  created  directly 
with  no  intermediate  steps.  The  flow  of  the  program  is  determined  by 
the  use  of  control  cards,  described  below. 

The  program  produces  graphs  very  similar  to  those  created  by  the 
Map  Display  program.  However,  the  graph  scale  used  by  the  Zonal  Pro¬ 
file  program  is  fixed,  so  that  the  profiles  of  variables  from  different 
experiments  can  be  compared  in  magnitude  as  well  as  in  shape.  The 
zonal  mean  values  are  printed  out  automatically  for  each  graph. 

PROFILE  DESCRIPTIONS 

The  zonal  arrays  written  by  the  M/A  program  are  used  to  form  32 
different  profile  graphs.  For  a  given  time  interval,  only  one  profile 
can  be  drawn  for  each  profile  number. 

The  formulas  describing  the  profiles  are  given  by  Table  12.  The 
mnemonics  used  are  given  in  Sec.  II,  and  represent  the  zonal  mean  ar¬ 
rays  after  they  have  been  multiplied  by  the  constants  given  in  Table  7. 
The  profile  titles  and  the  nominal  scale  factors  are  also  given  by 
Table  12. 

For  convenience,  some  of  the  profile  values  are  printed  in  two 
different  units.  For  profiles  1-5  and  7-16,  the  values  are  given  in 
terms  of  langleys  per  day  and  langleys  per  thousand  minutes  (KMIN  .  A 
sample  profile  graph  is  given  in  Fig.  3. 

INPUT  DESCRIPTIONS 

The  Zonal  Profile  program  uses  control  cards  input  by  the  user  to 
control  the  input  and  execution  of  the  program. 

Nine  types  of  control  cards  are  built  into  the  Zonal  Profile  pro¬ 
gram.  Only  the  first  four  columns  are  used  to  identify  the  control 
card,  with  the  other  76  columns  free  for  comments.  If  the  program  is 
expecting  a  control  card  and  does  not  recognize  the  next  card  as  a 
control  card,  the  card  is  still  printed,  and  three  lines  are  then 
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Table  12 

ZONAL  PROFILES  AVAILABLE 


Profile 

Scalea 

1 - - - 

Number 

Left 

Right 

Title 

Formula 

1 

0 

1600 

lnaolatlon  at  the  top  of  the  Atmoaphere 

(ly/day,  ly/kmln) 

scjlsz 

2 

0 

100 

Short-Wave  Abaorptlon  above  the  Tropoaphere 

(ly/day,  ly/kmln) 

STRA 

3 

0 

800 

Short-Wave  Abaorptlon  In  the  Tropoaphere 

(ly/day,  ly/kmln) 

AS1  +  AS  3 

4 

0 

800 

Total  Short-Wave  Abaorptlon  In  the  Atmosphere 

(ly/day,  ly/kmln) 

AS1  +  AS3  +  STRA 

5 

0 

1600 

Total  Reflection  In  the  Earth-Atmoaphere  System 

(ly/day,  ly/lcmln) 

RET  (IT 

6 

0 

1.0 

Albedo  (unltleaa) 

RFT0T/SC0SZ 

7 

0 

800 

Short-Wave  Abaorptlon  by  the  Earth 

(ly/day,  ly/kmln) 

S4 

8 

0 

800 

Effective  Long-Wave  Radiation  from  the  Earth 

(ly/day,  ly/kmln) 

R4 

9 

0 

800 

Long-Wave  Radiation  Loat  to  Space 

(ly/day,  ly/kmln) 

R0 

10 

0 

800 

Effective  Long-Wave  Radiation  Loat  frora  Atmosphere 

(ly/day,  ly/kmln) 

PC 

o 

1 

PC 

11 

-200 

600 

Heat  Loat  by  Evaporation  from  the  Earth 

(ly/day,  ly/kmln) 

E4  *  L/10 

12 

-400 

400 

Senalble  Heat  Loss  from  the  Earth 

(ly/day,  ly/kmln) 

F4 

13 

-300 

500 

Energy  Balance  of  the  Earth 

14 

(ly/day,  ly/kmln) 

S4  -  R4  -  E4  *  L/10  +  C1CE 

-800 

0 

Radiation  Balance  of  the  Atmoaphere 

(ly/day,  ly/kmln) 

R4  -  R0  +  AS1  +  AS3 

15 

-400 

400 

Energy  Balance  of  the  Atmoaphere 
(ly/day,  ly/kmln) 

R4  -  R0  +  AS1  +  AS  3  +  F4 

+  PPT0T  *  L/10  -  E4  *  L/10 

16 

-300 

500 

Energy  Balance  of  the  Earth-Atmoaphere  System 
(ly/day,  ly/kmln) 

S4  +  AS1  +  AS  3  -  R0  +  C1CE 
+  PPT0T  *  L/10 

17 

-1 

9 

Evaporation  (mm/day) 

E4 

18 

0 

10 

Precipitation  (mm/day) 

PPT0T 

19 

-5 

5 

Moisture  Balance  of  the  Atmosphere  (ram/day) 

E4  -  PPT0T 

20 

0 

4 

Latent  Heating  at  Sigma  3/4  (deg/day) 

HL3 

21 

0 

4 

Latent  Heating  at  Sigma  1/4  (deg/dsy) 

HL1 

22 

0 

4 

Short-Wave  Heating  at  Sigma  3/4  (deg/day! 

HSW3 

23 

0 

4 

Short-Wave  Heating  at  Sigma  1/4  (deg/day) 

HSW1 

24 

-4 

0 

Long-Wave  Heating  at  Sigma  3/4  (deg/day) 

HLW3 

25 

-4 

0 

Long-Wave  Heating  at  Sigma  1/4  (deg/day) 

HLW1 

26 

-2 

2 

Senalble  Heat  from  the  Earth  (deg/dsy) 

HF4 

27 

-2 

2 

Total  Heating  at  Sigma  3/4  (deg/day) 

HL3  +  HSW3  +  HLW3  +  HF4 
or 

HD1F  +  (HL3  +  HSW3  +  HLW3  +  HF4)/2.0 

28 

-  (HL1  +  HSW1  +  HLWl)/2.0 

-2 

2 

Total  Heating  at  Sigma  1/4  (deg/day) 

HL1  +  HSW1  +  HLW1 

or 

HDIF  +  (HL1  +  HSW1  +  HLW1J/2.0 

29 

0 

-  (HL3  +  HSW3  +  HLW3  +  HF4)/2.0 

1 

CL1  —  Deep  Convective  Cloudiness 

CL1 

30 

0 

i 

CL2  —  Large-Scale  Cloudiness 

CL2 

31 

0 

i 

CL3  —  Lose-Leval  Convective  Cloudiness 

CL3 

32 

0  . 

i 

CL  —  Total  Cloudiness 

CL 

PROFILE  19  EXPERIMENT  15  -  SAHARA  LAKE  -  SEPT  1972 

MOISTURE  BALANCE  OF  THE  ATMOSPHERE  ( MM/DAY  1  RRR 

TIME  =  370.00  TO  400.00  RANGE  IS  FROM  -5.00  TO  P.00 

NO.  fCAN  *  0.565  SO.  ft  AN  *  -0  508  GLOBAL  MEAN  *  0.057 

LAT  ZONAL  HANS 


90. 

-0.103 

86. 

-0.855 

85. 

-0.654 

78. 

-0.507 

74. 

-0.379 

70. 

-0.656 

66. 

-0.908 

65. 

-0.836 

58. 

-0.651 

54. 

-1.576 

50. 

-1.445 

46. 

-1.504 

45. 

-0.663 

38. 

0.195 

34. 

1.853 

30. 

3.008 

56. 

3.870 

55. 

3.504 

18. 

5.574 

14. 

1.669 

10. 

-0.151 

6. 

-0.316 

5. 

-1.585 

-5. 

-5.555 

-6. 

-5.638 

-10. 

-5.117 

-14. 

-1.763 

-18. 

-0.964 

-55. 

0.168 

-56. 

0.567 

-30. 

1.815 

-34. 

5.654 

-38. 

1.993 

-45. 

0.719 

-46. 

-0.455 

-50. 

-1.513 

-54. 

-1.659 

-58. 

-1.685 

-65 

-0.980 

-66. 

-0  479 

-70. 

-0.494 

-74. 

-0.569 

-78. 

-0. 154 

-85. 

-0.311 

-86. 

-0.410 

-90. 

-0.573 

Fig.  3 — Sample  graph  output  by  the  Zonal  Profile  program 
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skipped  in  the  printout  to  Indicate  the  possible  error.  This  continues 
until  a  control  card  is  found. 

The  control  card  DATA  controls  the  input  of  the  profile  descrip¬ 
tions.  For  each  graph  to  be  generated,  the  profile  number,  the  start 
and  end  times,  the  delta  time  (measurement  interval),  and  the  four- 
character  identification  code  word  must  be  input. 

To  simplify  the  DATA  card  input,  the  user  can  input  default  values, 
or  standard  values.  The  user  will  in  many  cases  wish  to  generate  sev¬ 
eral  profiles  for  the  same  time  period,  possibly  using  the  same  four- 
character  identification  code  word.  The  user  then  inputs  a  default  set 
of  time  values,  plus  a  default  code  word,  and  need  input  only  the  pro¬ 
file  number  for  the  profile  description.  The  code  word  serves  no  com¬ 
putational  purpose;  it  is  four  arbitrary  characters  that  the  user  can 
have  printed  on  each  graph  for  identification  purposes. 

The  delta  time  tells  how  often  to  plot  the  accumulated  profile 
values,  input  in  units  of  days,  so  that  more  L.,an  one  profile  can  be 
drawn  with  one  profile  description.  If  the  delta  time  is  zero,  then 
one  profile  is  drawn,  using  the  full  start  to  end  time  accumulation. 

If  snapshots  are  desired,  the  delta  time  must  be  equal  to  the  increment 
on  disk  (.25  day).  For  instance,  if  four  profiles  are  wanted  between 
370.0  and  371.0  days,  the  three  numbers  to  input  are  370.0,  371.0  and 


.25  days. 

A  running  count  of  the  number  of  profile  descriptions  is  kept,  so 
that  the  standard  values  can  be  changed  during  a  run.  Suppose  the 
standard  values  are  used  and  two  30-day  accumulations  are  desired. 

Then  the  standard  values  are  input  for  the  first  30  days,  and  the  DATA 
control  card  reads  in  the  descriptions  for  the  32  maps.  In  this  case, 
only  the  profile  numbers  need  to  be  input,  as  standard  times  and  code 
words  are  used.  Then  the  next  standard  times  for  the  next  30  days  are 
input  and  the  DATA  card  is  again  used,  again  with  only  the  profile 
numbers  1  through  32  on  the  cards.  This  procedure  has  input  the  data 
for  64  maps.  The  program  has  a  limit  of  128  profile  descriptions  for 
one  reading  of  the  M/A  disk,  but  these  sets  of  profiles  can  be  repeated 
as  many  times  as  desired  in  one  run. 


-36- 


Nine  control  cards  are  as  follows: 

TITLE :  The  card  following  this  has  the  run  title  in  columns  1-72. 
This  title  is  printed  on  each  subsequent  profile  and  printout.  Only 
one  title  is  used  for  each  set  of  profiles.  The  program  expects  a  con¬ 
trol  card  after  the  title. 

DATA :  This  card  tells  the  program  that  the  profile  descriptions 
follow.  The  profile  descriptions  are  terminated  by  a  blank  card,  after 
which  the  program  expects  another  control  card.  Each  of  the  five  data 
elements  uses  10  columns.  The  profile  number  is  an  integer,  right  ad¬ 
justed  in  columns  1-10.  The  start,  end,  and  delta  times  which  follow 
the  profile  number  are  floating-point  values.  The  code  word  is  input 
in  columns  47  through  50.  If  the  start  or  end  time  is  not  a  multiple 
of  .25  day,  or  the  delta  time  is  greater  than  .25  day  and  is  not  a 
multiple  of  .25  day,  an  error  message  is  written  and  the  profile  de¬ 
scription  is  deleted. 

The  number  of  input  profile  descriptions,  NMAP,  cannot  be  over  128, 
but  more  profile  descriptions  can  be  accumulated  again  after  the  128 
are  executed,  by  resetting  NMAP  to  zero  and  reading  in  more  data  cards. 

NMAP :  This  card  resets  the  number  of  input  profile  descriptions. 

The  number  is  input  right  adjusted  in  columns  1-10  of  the  next  card. 

If  a  second  set  of  profiles  is  to  be  generated  by  the  GO  control  card, 
the  number  of  descriptions  has  to  be  reset  to  zero.  If  the  NMAP  card 
is  used  to  input  a  negative  number,  then  the  number  of  profile  descrip¬ 
tions  is  reduced  by  that  amount.  The  program  expects  a  control  card 
after  the  one  data  card. 

CODE:  This  card  sets  a  default  four-character  code  word  for  the 
DATA  cards.  If  the  code  word  for  the  data  description  is  blank,  this 
word  is  inserted.  It  will  usually  be  used  to  input  the  initials  of 
the  interested  researcher  for  identification  purposes.  The  code  word 
is  input  in  columns  1-4  of  the  next  card. 

STANDARD :  The  default  start,  end,  and  delta  times  are  read  in 
following  this  card.  The  data  elements  use  10  columns  each,  and  are 
floating-point  values,  starting  in  column  1. 

SCALES:  This  control  card  changes  at  least  one  of  the  default 
profile  scales.  For  a  specific  experiment,  the  built-in  scales  of  the 
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graph  area  may  not  be  appropriate,  and  they  can  be  changed  here.  One 
new  scale  is  input  per  card;  the  data  cards  are  followed  by  a  blank  card. 
The  profile  type  number  is  input  as  a  right-adjusted  integer  in  columns 
1-10,  the  left  graph  value  in  columns  11-20,  and  the  right  graph  value 
is  columns  21-30.  Both  values  are  input  as  floating-point  numbers. 

GO:  This  control  card  tells  the  program  to  execute  the  profiles 
that  have  been  input.  No  data  cards  are  input  for  this  control  card. 

The  profile  descriptions  are  printed  out,  and  then  the  graphs  and 
printout  generated.  After  all  the  profiles  are  done,  the  program  re¬ 
turns  to  process  more  control  cards. 

NUMZON :  This  control  card  determines  the  number  of  zonal  mean  ar¬ 
rays  in  the  input  data  sets.  Default  is  24  arrays.  Note:  Experiments 
10-13  have  24  arrays,  experiments  14-17  have  25  arrays. 

STOP:  The  program  must  end  with  this  control  card  for  normal 
termination.  No  cards  are  read  beyond  this  control  card. 

Table  13  gives  a  short  description  of  each  of  the  control  cards, 
with  the  number  and  format (s)  of  the  data  cards  which  follow  each. 


Table  13 

CONTROL  WORDS  FOR  THE  ZONAL  PROFILE  PROGRAM 


Control 

Word 

Description 

Number  of 
Data  Cards 

Format  of 
Data 

TITLE 

Run  title 

1 

18A4 

DATA 

Profile  descriptions — No.,  times,  ID 

N 

110,  6X,  A4 
3F10.0 

C0DE 

Standard  4-character  identification 
code  word 

1 

A4 

STANDARD 

Standard  start,  end,  and  delta  times 

1 

3F10.0 

SCALES 

Change  profile  scales 

N 

110,  2F10.0 

Print  out  input  and  execute 

0 

None 

NMAP 

Change  the  number  of  input  profile 
descriptions 

1 

110 

NUMZ0N 

Set  the  number  of  zonal  means 

1 

10 

ST0P 

End  program  normally 

0 

None 
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SAMPLE  INPUT  RUN 

The  card  input  listing  in  Fig.  4  is  a  valid  input  deck  for  the 
Zonal  Profile  program.  First  the  title  is  input,  followed  by  the  stan¬ 
dard  times  and  code  words.  The  data  will  be  accumulated  from  day  340 
to  day  370,  arbitrarily  using  the  identification  code  STD,  standing 
for  the  standard  output.  The  DATA  control  card  then  inputs  all  32  pro¬ 
files,  using  the  standard  values.  Note  that  only  the  profile  numbers 
are  input  in  this  section. 

As  the  rest  of  the  maps  for  this  pass  are  not  standard,  the  code 
word  is  changed  from  STD  to  DSP.  Then  three  profile  descriptions  are 
input;  thus  35  profile  descriptions  have  been  input.  The  first  of  the 
new  descriptions  draws  five  graphs  of  profile  No.  5  at  the  times  399.0, 
399.25,  399.5,  399.75,  and  400.0.  As  the  M/A  program  accumulates  these 
values,  the  maps  represent  six-hour  accumulations,  even  though  only  one 
time  point  is  read  from  disk  for  each. 

The  next  profile  description  generates  one  graph,  and  accumulates 
the  data  for  profile  8  between  395  and  400  days,  with  the  code  word 
WLG  used.  Again,  the  accumulation  averages  the  dick  time  points  395.25 
to  400.0. 

The  last  profile  description  generates  nine  graphs  of  profile  10, 
using  the  nonstandard  code  word  DSP.  The  times  are  390  to  391,  etc. 

Each  accumulation  uses  four  disk  time  sets  for  each  dav. 

The  control  card  GO  then  outputs  the  input  profile  descriptions, 
the  scales,  etc.,  and  the  profiles  are  generated.  After  the  profiles 
are  drawn,  control  is  again  returned  to  the  input  cards.  The  number 
of  profile  descriptions  is  reset  to  zero,  and  new  standard  values  are 
input.  In  this  case,  the  values  are  to  be  accumulated  for  30  days, 
between  day  400  and  day  430,  using  the  code  name  STD  for  all  of  the  32 
profiles.  These  profiles  are  then  executed  with  the  GO  control  card, 
and  the  program  ends  with  the  STOP  control  card. 
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VI.  USER’S  GUIDE  FOR  THE  QUICK  LOOK  PROGRA*  1 


Quick  Look  is  designed  to  help  the  programmer  and  researcher  to  de¬ 
cide  if  the  M/A  program  is  running  correctly.  At  checkpoints  through¬ 
out  the  simulation,  the  M/A  program  may  be  stopped  and  the  data  on  disk 
irspected  to  determine  if  the  program  is  behaving  as  expected.  The 
Qrick  Look  program  was  designed  primarily  to  do  this. 

Quick  Look  reads  both  the  large  array  data  set  and  the  zonal  mean 
data  set.  The  program  also  has  a  limited  capability  of  accumulating 
values’,  the  three  precipitation  arrays  are  accumulated  to  form  the 
total  precipitation  values,  and  all  the  zonal  mean  arrays  can  be  accum¬ 
ulated.  Of  course,  logic  can  be  added  to  accumulate  any  of  the  arrays 
on  the  disk,  but  only  those  mentioned  are  automatically  implemented. 

At  a  specified  time  point,  the  program  will  print  out  all  the  re¬ 
start  arrays:  P,  Ul,  U3,  VI,  V3,  Tl,  T3,  Q3,  GW,  and  GT.  The  user 
inputs  the  start  and  end  times  of  the  disk,  and  at  every  Nth  time  point, 
where  N  is  input  by  the  user,  values  at  a  specified  grid  point  from 
each  large  array  are  printed  together  with  all  the  zonal  mean  arrays 
(accumulated  from  the  previous  six  hours). 

CARD  INPUT 

The  first  card  of  the  input  is  a  title,  80  columns  long.  The 
next  cards  describe  the  printout.  The  start  and  end  times  are  piven 
in  columns  11-12  and  13-24,  respectively.  Time  is  input  in  days,  and 
must  be  exactly  the  same  as  on  the  disk.  Columns  25-36,  right  adjusted, 
give  the  integer  number  N,  which  controls  the  printout  of  the  large 
arrays  at  the  specified  grid  point  and  the  snapshot  of  the  zonal  mean 
arrays.  Columns  37-48,  right  adjusted,  hold  the  integer  flag  for  the 
accumulations  over  the  entire  time  period.  If  it  is  zero,  no  accumu¬ 
lations  are  done;  if  not  zero,  the  zonal  means  are  accumulated  and 
printed  out.  If  the  integer  flag  value  is  above  9,  then  the  total 
precipitation  values  are  also  accumulated. 
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The  time  at  which  to  print  all  the  restart  arrays  is  set  within 
the  program  by  a  variable  (PRTALL) .  This  must  be  set  by  the  programmer 
or  researcher  before  compilation. 

The  Quick  Look  program  executes  the  one-card  instruction,  and  then 
looks  for  another  data  card.  The  one  limitation  on  this  following  card 
is  that  its  starting  time  must  be  greater  than  the  ending  time  on  the 
previous  card.  The  program  is  terminated  by  a  blank  card. 

PRINTOUT  DESCRIPTION 

Suppose  the  start  and  end  times  are  405  and  410,  N  is  set  to  200, 
and  the  flag  set  to  10.  Then  at  day  405,  the  array  grid-point  values 
are  printed  out,  followed  by  the  zonal  mean  values  accumulated  for  the 
previous  six  hours.  The  grid-point  values  are  printed  out  at  -14°  lat¬ 
itude  and  35°  W  longitude,  corresponding  to  the  latitude  grid  number 
20  and  longitude  grid  number  30.  (This  grid  point  can,  of  course,  be 
easily  changed  within  the  program  to  any  other  grid  point.)  Since  N 
was  set  to  200,  this  printout  would  be  repeated  after  200  disk  time 
points,  or  50  days,  but  as  the  end  time  come  before  then,  this  printout 
is  done  only  at  the  start  time,  day  405. 

Supposing  that  the  print  time,  PRTALL,  was  set  to  410,  all  the 
restart  arrays  and  precipitation  arrays  are  printed  out  at  that  time. 
Each  array  takes  eight  pages  to  print  out;  an  eight-character  identifi¬ 
cation  for  each  array  is  built  into  the  program  and  is  printed  at  thu 
top  of  each  page.  Next  the  accumulated  zonal  arrays  are  printed  out, 
preceded  by  the  title  and  the  start  and  end  times  of  the  accumulation. 
The  24  zonal  arrays  are  printed,  each  value  preceded  by  its  latitude, 
with  a  built-in  four-character  identification  given  as  a  heading. 

After  this,  the  precipitation  array  is  printed.  Note  that  the  zonal 
array  values  are  always  printed  out  after  being  multiplied  by  the  con¬ 
version  constants  given  in  Table  7. 

Two  subroutines  are  included  in  the  Quick  Look  program  to  do  the 
printing  of  the  large  arrays  and  the  zonal  means.  The  array  and  a 
snort  name  are  input  to  the  routines  and  the  data  are  printed  out. 

The  printout  of  the  large  arrays  gives  the  latitude  and  longitude  of 
the  rows  and  columns.  Each  data  value  is  printed  with  five-place  ac¬ 
curacy.  The  zonal  arrays  are  printed  next  to  the  latitude  values. 
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As  a  bonus,  other  data  are  automatically  printed  out.  The  first 
time  that  TAU  and  the  constant  C-array  (see  Sec.  II)  are  read,  all  800 
constants  of  this  C-array  are  printed  out  in  three  forms:  as  integers, 
binary-coded  decimal  characters,  and  floating-point  values.  The  24  or 
25  constants  used  in  converting  zonal  mean  arrays  are  also  printed  out. 

By  adding  one  card,  a  feature  can  be  used  which  reallv  helps  the 
researcher  and/or  programmer  analyze  the  output.  The  topography  array 
T0P0G  gives  the  sea,  land,  and  ice  areas  of  the  Earth.  Within  the  pro¬ 
gram,  this  array  is  changed  by  setting  0.0  to  represent  ocean,  1.1111 
to  represent  ice,  and  9.9999  to  represent  land.  The  resultant  array 
can  be  printed  by  adding  one  card,  so  that  the  user  can  determine  at  a 
glance  whether  a  particular  grid  point  is  ocean,  land,  or  ice. 

JOB  CONTROL  COMMENTS 

The  Quick  Look  program  uses  less  than  150K  bytes  of  core,  and  takes 
less  than  two  minutes  to  run  on  the  360/91.  The  number  of  pages  gen¬ 
erated  may  be  very  large,  depending  on  what  is  being  printed  out.  The 
restart  disk  data  set  is  input  on  unit  10,  and  the  zonal  disk  data  set 
is  input  on  unit  11.  By  concatenating  df.ta  sets,  20  days  or  more  can 
be  investigated  on  one  pass. 
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VII.  USER'S  GUIDE  FOR  THE  DISK  INITIALIZER  PROGRAM 

The  Disk  Initializer  program  is  used  in  conjunction  with  the  M/A 
program.  It  will  read  a  time  history  up  to  a  specified  time  point, 
then  write  the  data  onto  a  new  disk  so  that  the  M/A  program  can  continue 
its  simulation.  Because  of  the  way  the  M/A  program  is  written,  it  is 
easier  and  safer  to  use  this  program  rather  than  the  M/A  program  itself 
to  switch  disks  when  writing  the  time  history  data  sets. 

Disk  Initializer  has  other  uses,  too.  It  can  easily  be  modified 
to  set  up  the  initial  conditions  of  a  new  experiment.  Suppose  the  new 
experiment  modifies  the  topography  array;  then  the  Initializer  can  read 
in  the  new  values  and  initialize  the  new  disk  using  these  values.  For 
printout  purposes,  the  subroutine  in  the  Quick  Look  program  that  prints 
out  the  large  arrays  can  be  copied  and  added  to  Disk  Initializer  to 
print  out  the  old  and  new  topography  arrays. 

Disk  Initializer  is  useful  in  emergency  situations.  No  computer 
is  perfect,  and  when  a  program  uses  30  hours  of  CPU  time,  it  is  exposed 
to  a  chance  of  malfunction.  Sometimes  a  special  emergency  program  has 
to  be  written  to  save  all  or  part  of  an  M/A  disk;  Disk  Initializer  can 
be  used  as  a  basis  for  the  new  program. 

Disk  Initializer  uses  400K  of  core  and  less  than  one  minute  of  CPU 
time.  The  old  restart  data  set  is  on  unit  11,  the  old  zonal  data  set 
is  on  unit  13,  the  new  restart  data  set  is  on  unit  12,  the  new  zonal 
data  set  is  on  unit  14,  and  the  new  last-time  record  is  on  unit  1. 

The  time,  in  days,  to  read  up  to  and  write  on  the  new  units  is  set  by 
the  variable  DAYS  within  the  program  itself.  There  are  no  cards  input 
to  this  program. 
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VIII.  USER'S  GUIDE  FOR  THE  DIFFERENCE  PROGRAM 


The  Difference  program  is  used  in  conjunction  with  MAPGEN  and  Map 
Display  for  the  purpose  of  comparing  two  maps.  Unlike  the  other  two 
auxiliary  programs,  the  Difference  program  is  relatively  stable  and 
will  not  normally  change  from  experiment  to  experiment. 

Two  tapes  from  MAPGEN  runs  -re  used  as  input  to  generate  a  thi rd 
tape  giving  the  differences  of  values  on  the  original  two  tapes.  This 
third  tape  can  then  be  used  as  input  to  the  Map  Display  program. 

The  card  input  is  identical  to  the  MAPGEN  card  input,  in  format: 
a  title  card,  followed  by  data  sets  consisting  of  the  first-card  time 
description,  followed  by  the  map  descriptions.  The  map  descriptions 
do  not  have  to  be  in  the  same  order  as  in  MAPGEN,  and  the  limit  of  40 
map  descriptions  per  set  does  not  apply. 

The  values  differenced  are  the  large  46  *  72  data  arrays,  the 
zonal  means,  the  zonal  means  over  the  ocean,  and  the  global  means.  The 
values  from  the  "low"  tape,  the  one  on  FORTRAN  unit  10,  are  subtracted 
from  the  values  from  the  "high"  tape,  on  unit  9,  and  the  differences 
put  on  the  output  tape  on  unit  11.  Only  the  map  descriptions  input  on 
cards  are  written  on  the  tape.  To  match  the  two  sets  of  values,  the 
start  and  end  times,  the  map  number,  and  the  SURF  values  must  be  the 
same.  For  each  map,  the  values  written  on  the  output  tape,  such  as  the 
ID,  the  map  title,  and  the  number  of  accumulations,  are  the  same  as  on 
the  high  tape. 

The  Difference  program  uses  less  than  150K  bytes  of  core  and  takes 
less  than  two  minutes  of  CPU  time  on  the  IBM  360/91. 
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APPENDIX  A 

LISTING  OF  STANDARD  INPUT  PACKAGE 
FOR 

MAPGEN  AND  MAP  DISPLAY 


THIS  IS  THE  ST ANOARO  INPUT  PACKAGE  FOR  NAPGEN  AND 
NAP  0 1  SPLAY •  COLUMNS  1-7  -  NAP  NUMBER  COLUMNS  3-7  - 
SURF  COLUMNS  B-l?  -  CONTOUR  INTERVAL  COLUMNS  13-17  - 
OOTTFD  LINE  COLUMNS  73-77  -  NO.  OF  PRINTOUTS  COLUMNS 
2B-3?  _  FLAG  FOR  USERS  CONTOUR  VALUES  COLUMNS  33-37  - 
TONAL  MAP  FLAG  COLUMNS  3B-A?  -  CONTOUR  MAP  FLAG 

THE  FOLLOWING  FORMS  ALL  3  TVPES  OF  OUTPUT  IN  MAP 
DISPLAY. 
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Appendix  C 

CORRESPONDING  MAP  NUMBERS  IN  THE  OLD3  AND 
NEW  VERSIONS  OF  MAPGEN 


"Old”  map  numbers  are  those  used  in  the 
earlier  documentation  1). 


-56- 


Table  D1 

SYMBOLS  USED  IN  MAPGEN  MAIN  ROUTINE 

Name  Meaning 

AVE  Contour  Interval  for  Map  Display 
AVFR(46 , 72 ,40)  Accumulation  arrays  for  current  data  set 

C0N0  Flag  for  Input  contour  intervals  for  Map  Display 
C0NST  Value  of  weighting  factor  for  this  time  point 

Dl.N  Dashed-llne  lsollne  value  for  Map  Display 

EACH  No  longer  used 

EXPN  Flag  for  punched  output,  set  to  the  Integer  NEXP 
HE\D(19)  Run  title 

ID  Temporary  storage  for  4-character  Identification 

word 

ISRC  Flag  for  M/A  Input  data  type--set  to  new  disk 
format 

KEYS(32)  Flags  not  used  In  MAPGEN 

LEV  Flag,  *  true  If  SURF  is  a  sigma  level 

LHLF  Flag  used  In  some  map  routines;  «  true  If  SURF 

greater  than  0.5 

MAP(41)  Map  numbers  for  current  data  set 

NAME(13)  Temporary  storage  for  map  titles 

NEXP  Integer  value  of  EXPN 

N00UT  Flag  for  C0MP3;  call  *  0  if  C0MP3  not  called  for 
the  time  point 

PRN0  Print  flag  for  Map  Display 

SNT(40)  Storage  for  variables  no  longer  uspd  In  MAPGEN 

SRF(40)  SURF  values  for  current  data  set 

SURF  The  sigma  level,  pressure,  or  flag  describing  the 

map 

TEND  Ending  time  of  accumulation  of  data  set 

TICE  The  melting  point  of  ice  (-  273.1  deg  K) 

TO  Starting  time  of  accumulation  of  data  set 

WF(46,2)  Weighting  functions  for  zonal  means 

WORK2(46,92)  The  map  values  for  the  current  time  point 

XID(40)  Storage  foi  4-character  Identification  word  for  maps 

XNAMF(13,40)  Storage  for  map  titles 

XN0(4O)  Storage  of  weighting  factors  for  accumulations 

XSTAGI(40)  Flags  for  r  or  u,v  grids  In  zonal  direction 

XSTACJ(40)  Flags  for  »  or  u,v  grids  In  meridional  direction 

XZM(46,40)  Accumulation  of  zonal  means 

XZMM(40)  Accumulation  of  global  means 

X7V2(46,40)  Accumulation  of  zonal  means  over  oceans 

ZM(46)  Temporary  storage  for  zonal  means 

ZMM  Temporary  storage  for  global  mean 

ZM2(46)  Temporary  storage  for  zonal  means  over  oceans 

ZN0  Value  of  weighting  factor  for  current  map 

accumulation 
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Table  D2 

THE  LARGE  MAPGF.N  ARRAYS 


Name 

Dimensions 

Description,  Use,  and/or  Map  Reference 

A1 so  Cal  led 

Stored  in 

ACPPC1 

(46,72) 

Accumulated  precipitation  arr.iv--soe  user's  Guide 

.... 

Ml  SC 

ACPPC3 

(46,72) 

Accumulated  precipitation  array — see  User’s  Guide 

—  - 

M1SC 

ACPPLS 

(4b, 72) 

Accumulated  precipitation  array-  see  User’s  Guide 

— 

M1SC 

ACRET0 

(46,72) 

Sum  of  RET0T  arrays  —  see  Map  28 

— 

DATA28 

ACSC0S 

(46,72) 

Sum  of  SC0SZ  arrays — see  Map  28 

— 

DATA28 

ALS  1  n 

ALBED0 

(46,72) 

Surface  albedo 

C0MP3  and 
RET0T  cal- 

Mi  SC 

cul at  ions 

AVER 

(46,72,40) 

Sum  of  map  data  arrays 

— 

Local  variable 
in  MAiN 

C  C 

(46,72,4) 

Cloudiness:  CL1  ,  CL2 ,  CL3,  and  CL— see  Map  14 

... 

EXC0M 

C0NM 

(46,72) 

Temporary  storage  for  Map  7,  not  In  any  common 

... 

Local  variable 
MAP  7  rout i nr 

C0NV 

(46,72) 

Packed  S4  and  T4 — see  Maps  9  and  21 

SD 

Ml  SC 

CPC1 

(46,72) 

Latent  heating  in  layer  1,  Cl  +  PCI— see  Map  24 

... 

EXC0M 

CPC  3 

(46,72) 

Latent  heating  in  layer  3,  C3  +  P3 — see  Map  24 

I.XC0M 

FD 

(46,72) 

Temporary  storage.  Map  7  and  C0MP  3 

PL 

H(  .6,72,.  ) 

MISl 

TD 

GT 

(46,72) 

Ground  temperature,  part  of  the  state  vector 

... 

•JAR  V 

CW 

(46,72) 

Ground  wetness,  part  of  state  vector 

— 

•JARS 

(46,72,1) 

Temporary  storage  for  Map  7 

1  PHI 

PV 

H 

|w 

MISC 

(46,72,2) 

Temporary  storage  Map  7  and  C0MP3 

PC 

FD 

TD 

P 

(46,72) 

Surface  pressure  parameter,  part  of  state  vector 

— 

QARV 

PHI 

(46,72,1) 

Temporary  storage  for  Map  7 

ri(s6,72.l) 

l*V 

W 

MISC 

PRCLH 

(46,72) 

Latent  heating  in  layers,  ('RFC  •  C|.H — see  Map  24 

... 

EXCOM 

PREGT 

(46,72) 

Pre  ‘'MPT  ground  temperature 

... 

PRE  3 

PREGW 

(46,72) 

Pre-C0MP3  ground  wetness 

... 

PRE  1 

PRFQ3 

(46,72) 

Pro-C0MP3  lovol  \  mixing  ml  in 

... 

PRE  1 

I'RKT 

(46,72,2) 

Pn-CflMPi  lev**  Is  1  and  1  temporal ures 

... 

PRE  ) 
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Naae  Dimensions  Description,  Use,  and/or  Map  Reference 


Also  Called!  Stored  In 


PT  (46,72)  Packed,  HI  and  H3,  total  heating — see  Map  27 


PU  (46,72)  Temporary  storage,  Map  7  and  C0MP3 


PV  (46,72)  Temporary  storage  for  Map  7 


I  H(46 ,72 ,2)  M1SC 

[td  1 

H(46  ,72 , 1) 

PHI  MISC 

W 


Q3  (46,72)  Level  3  mixing  ratio,  part  of  state  vector 

Q3T  (46,72)  CICE,  heat  rising  through  Ice — see  Map  23 

RET0T  (46,72)  Total  reflected  solar  radiation — see  Map  28 

SD  (46,72)  Packed  S4  and  T4 — see  Maps  9  and  21 

SN  (46,72)  Packed  R2  -  R0,  R4  -  R2— see  Mbp  26 


VT(46 ,72 ,2)  QTARY 


SSC0SZ  (46,72) 


Solar  radiation  at  top  of  atmosphere, 
SC0SZ — see  Map  28 


(46,72,2)  Levels  1  and  3  temperatures,  part  of  state  vector 


TD  (46,72)  Temporary  storage,  Map  7  and  C0MP3 


H(46,72,2)  MISC 
PU 


T0P0C  (46,72) 


TS  (46,72,2) 


[('*6,72,1) 


(46,72,2) 


Topography:  temperature  positive  for  ocean,  less 
than  10"®  for  land  Ice 

R4  net  upward  long-wave  radiation  at  surface — 
see  Map  22 

Packed  AS1,  AS3,  solar  absorption  In  layers — 
see  Map  25 

Packed  F4,  E4,  evaporation  aad  sensible  heat — 
see  Maps  18-20 


UT(46 , 72  ,2)  QTARY 


U  (46,72,2)  Zonal  winds  for  levels  1  and  3,  part  of 


state  vector 


(46.72.1)  Temporary  storage  In  AVR7C 

(46.72.2)  R4  net  upward  long-wave  radiation  at  surface — 

see  Hap  22 

V  (46,72,2)  Meridional  winds  for  levels  l  and  3,  part  of 


state  vector 


(46.72.1) 

(46.72.2) 


Temporary  storage  in  map  subroutines 
Packed  R2  -  R0,  R4  -  R2— see  Map  26 


U  (46,72)  Temporary  storage  for  Map  7 

W0RK1  (46,72)  Temporary  storage  In  AVRX 

W0RK2  (46,72)  Temporary  storage  In  map  sul 


H(46 ,72 ,1; 

PHI  MISC 

PV 

UT(46 ,72 ,1)  QTARY 
VT(46 , 72 ,1)  QTARY 
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Table  D3 

MAP  DISPLAY  VARIABLES  IN  MAIN  ROUTINE 


Name 

Description 

AVE 

Input  contour  interval,  SINT  set  to  it 

AVO 

Dashed-line  value,  set  from  DLN 

CN 

Current  contour  Interval,  initialized  at  the  lowest 
contour  value 

CN2 

The  smallest  grid  value 

CN3 

The  largest  grid  value 

CONFLG 

Flag  for  drawing  the  contour  map 

CONIN(15) 

The  input  contour  values,  if  any 

CONINT(3) 

BCD  information  -  ' ISOLINES  AT' 

CONO 

Input,  real  number  of  contour  intervals;  zero  if  none 

CX(1500) 

x-axis  coordinates  of  lines  forming  the  continental 
outlines — each  set  of  two  points  defining  a  line 

CY (1500) 

y-axia  coordinates  of  lines  forming  the  continental 
outlines— each  set  of  two  points  defining  a  line 

DLAT 

Number  of  degrees  per  grid  point  on  y-axis 

DLN 

Input  dashed-line  value,  AVO  set  to  it 

DLON 

Number  of  degrees  per  grid  point  on  x-axis 

DT 

Used  in  an  older  version  of  Map  Display;  should  be 
input  as  zero  in  this  version.  At  one  time,  MAPGEN 
did  not  accumulate  the  arrays,  but  put  each  array 
on  tape  (one  for  each  map  every  6  hours,  for  usual 
case).  Map  Display  then  accumulated  the  arrays. 

The  option  was  deleted,  but  some  of  the  logic  still 
exists  in  Map  Display. 

DTT 

Small  time  increment  used  in  time  comparisons 

DX 

Number  of  rasters  corresponding  to  one  grid  interval 
on  the  x-axls  In  the  contour  map 

DY 

Nianber  of  rasters  corresponding  to  one  grid  interval 
on  the  y-axls  in  the  contour  map 

EACH 

Not  used;  should  be  input  as  zero 

EXPN 

Not  used 

FLAG  7 

Logical  flag  for  the  meridional  mass  flux  profile 
map  done  in  map  7 

FREQ 

Not  used,  but  either  this  or  EACH  should  be  input  as 
zero.  Input  slong  with  TA  in  columns  41-50 

IC 

Nun  be r  of  points  fa  CX,  CY  array.  IC/2  lines  drawn 

IM 

Integer  number  of  grid  points  on  x-axls,  set  to  XMAX 

JM 

Integer  number  of  grid  points  on  y-axls,  set  to  YMAX 

LAT(47  ) 

Latitude  grid  points  in  degrees,  used  in  POLP 

LIN(3) 

Used  in  CONTUR  to  speed  up  the  ieolins  drewing  in 

ICS.  In  COmON/SPEED/ 

MAP 

Map  number  read  from  cards 
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Table  D3— continued 


Name 


MAPNO 
NAME  7 

NOPR 

NOR 

pr 

POI(2) 

PRNO 

RAD 

REW 

SINT 

SINTM 

SRF 

SURF 

TA,  TD 

TB,  TK 
TC 

TEND 

TOT 

TO 

T1 

T2 

WF(46,2) 

WORK1 (46, 72) 

WORK2  (46,  72) 

WRK (46,  72) 

XMAX 

XI 

Y1 

ZMM 

ZM(46) 

ZM2(46) 

7.M3  (46) 

ZNCR 

ZNT  (3) 

7.0  PI. 


Description 


Map  number  read  from  tape 

ThLpCmiT  f°r  thC  merldIonal  mass  Hux  Profile 

Integer  number  of  times  to  print,  set  from  PRNO 
Integer  flag  for  projection  type;  set  from  ZNOR 

degrees^ '  ^  ^  convcrslon  of  radians  and 

No  longer  used 

Input,  real  number  of  times  to  print  the  map  data 
"  6‘375  '  10  meters ,  the  radius  of  the  earth 

Rewinefsear^°iSt°f  ^  'T  beln*  rewound  twlce 
while  searching  for  a  map’s  data  record 

Contour  interval,  set  from  AVE 

100  times  SINT,  the  contour  interval;  used  to  prevent 
drawing  of  more  than  100  contour  intervals 
SURF  value  read  from  card 
SURF  value  read  from  tape 
Starting  time  of  map  description,  In  days 
Ending  time  of  map  description,  in  days 
A  large  time,  set  to  100  days  beyond  TB 
End  time  read  from  tape,  in  hours 

*'u  n~«>*  PM...  , 

Start  time  read  from  tape,  in  hours 
Start  time  read  from  tape,  in  days 
End  time  read  from  tape,  in  days 

Input  continental  outline  data 
Urge  array  input  from  tape,  equivalenced  to  WRK 
rge  array  input  from  tape,  equivalenced  to  W0RK2 
Number  of  grid  points  on  x-axis,  floating  point 
Start  of  the  plotting  surface  on  the  x-axis 
Start  of  the  plotting  surface  on  the  y-axis 
Global  mean,  recalculated  from  tape  value 
Zonal  means  for  large  array,  recalculated 
Zonal  means  of  the  oceans  from  tape 

ArInym^7he  mCr‘dl0nal  maSS  flux  calculated 

°f  pr°Jectlon  contour 
nag,  NOR,  its  integer  value,  is  used 

BCD  information  -  ’INTERVAL-’ 

Flag  for  drawing  the  zonal  means  map 
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Table  D4 

ZONAL  PROFILE  PROGRAM  VARIABLES 
Name  Description 

ACCZON(46, 128)  Storage  for  time  accumulation  of  profiles 

CONSTS  In  ACCVAL,  the  current  weighting  factor  for  the  map 
data  added  to  ACCZON 

CONV(30)  Conversion  factors  for  the  Input  zonal  arrays 

FLGMAP  Flag  Indicating  that  a  map  is  to  be  plotted,  in  EXEQ 

ICARD(20)  Used  to  read  the  control  cards 
IPRO(129)  Profile  type  number  for  maps 
ITITLE  Title  for  maps 

NAMCOD(129)  Four-chavacter  identification  for  each  map 
NAMTYP(9)  Array  containing  the  control  words 

NCODE  Default  value  for  NAMCOD.  Normally  blank 
NMAP  Number  of  profile  maps  input — counter 
NPTZON  Number  of  zonal  divisions,  nominally  46 
NTYPES  Number  of  different  control  words 

NUM(129)  Number  of  accumulations  in  ACCZON  array 
NUMUNT  Unit  number  of  input  zonal  means 

NUMMAP  Number  of  types  of  profile  maps 

NUMZON  Number  of  input  zonal  means 

SCALE(2 , 32)  Range  for  profile  type  in  map  area 

SUMWT  The  weighting  factor  for  the  array  ACCZON,  in  EXEQ 

T(3,129)  Start  time,  end  time,  and  time  range  per  profile 
(input  in  days,  changed  to  hours) 

TAU  In  EXEQ,  time  of  disk  data  sets  in  hours 
TAULST  Time  of  last  accumulated  map  data 
TIMDIF  Time  increment  of  data  sets  on  disk 
TIMMAX  Maximum  time  demanded  in  profile  descriptions 

TIMMIN  Minimum  time  demanded  in  profile  descriptions 

TISTD  Default  value  for  start  time.  Nominally  0.1 

T2STD  Default  value  for  end  time.  Nominally  0.1 

T3STD  Default  value  for  time  range.  Nominally  0.4 

WT(46)  Weighting  function  used  in  calculating  zonal  means 

Z(200)  Mode  set  array  for  graphics  programs 
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Table  D5 

SUBROUTINE  MAPZON  VARIABLES 


Name 

Description 

DEL 

Number  of  rasters  from  ZPT  for  the  latitude  labels 

DELTA 

Number  of  rasters  to  move  label  If  two  columns  of  num¬ 
bers  are  printed 

DIFFG 

Number  of  rasters  in  the  plotting  surface  in  the  x  di¬ 
rection 

EDG 

Length  of  the  lines  bracketing  the  plotting  surface 

FLCINT 

Flag,  ■  true  if  the  accumulation  is  more  than  one  time 
point 

FLGMIN 

Flag,  -  true  if  two  columns  of  numbers  are  printed  in 
the  profile  map  (the  langleys/minutes  column) 

GRIDX1 

Left-hand  x-axis  value  of  the  plotting  surface 

GRIDX2 

Right-hand  x-axls  value  of  the  plotting  surface 

GRIDY1 

Lower  y-axia  value  of  the  plotting  surface 

GRIDY2 

Upper  y-axis  value  of  the  plotting  surface 

IDUM 

Duniny  variable  used  in  call  to  LINSEG 

IMAP 

Profile  map  number  to  plot 

IN0 (3) 

BCD  information,  'NO.  MEAN  »' 

IS0(3) 

BCD  information,  'SO.  MEAN  •' 

IZ0N(3) 

BCD  information,  'ZONAL  MEANS' 

MAPNO 

Proflle  type  number 

NTIT (17) 

Title  for  profile  type 

PLATP 

Latitude  counter  in  loop 

SCALE1 

Values  at  the  left  edge  of  the  plotting  surface 

SCALE2 

Values  at  the  right  edge  of  the  plotting  surface 

SUMWT 

Weighting  factor  for  ACCZON 

TAU 

Last  time  point  accumulated  for  this  profile 

T1 

Start  time  oi  accumulation  for  this  profile 

T2 

End  time  of  accumulation  for  this  profile 

VAL 

Raster  point  along  the  x-axls  to  print  latitude  labels 

VALCBL 

The  global  mean  for  the  profile  values 

VALN0 

Northern  hemisphere's  zonal  mean 

VALS0 

Southern  hemisphere's  zonal  mean 

XI, Y1 

Starting  point  of  zonal  mean  line  within  plotting  surface 

X2,Y2 

End  point  of  zonal  mean  line  within  plotting  surface 

YLN1 ,  to  5 

Raster  values  of  the  y-axls  for  lines  1  (top  line)  to 

5  of  the  headings 

YPT 

Raster  value  of  y-axis  in  drawing  graph 

ZPT 

Raster  value  of  zero  in  the  x-axis  in  the  plotting  sur¬ 
face 

Z2ARY (46) 

Zonal  mean  values  in  units  of  ly/min,  if  used 
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